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Abstract. Following work of Dupire (2005), Carr-Lee (2010) and Cox-Wang (2011) on connections between 
Root's solution of the Skorokhod embedding problem, free boundary PDEs and model-independent bounds on 
options on variance we propose an approach with viscosity solutions. Besides extending the previous results, it 
gives a link with reflected FBSDEs as introduced by El Karoui et. al. (1997) and allows for easy convergence 
proofs of numerical schemes via the Barles-Souganidis (1991) method. 



1. Introduction 

Let B be standard, real-valued Brownian motion (started at 0) and fx a centered integrable probability 
measure on (R,£>(R)). In 1961 Skorokhod [36] asked whether it is possible to find a stopping time r such that 
B T ~ fx and B T = (B TAt ) t is uniformly integrable; a stopping time r with these properties is now usually referred 
to as a solution of the Skorokhod embedding problem or Skorokhod stopping problem. Further, he provided an 
affirmative answer by constructing a stopping time which depends on an additional random variable independent 
of B. Skorokhods' motivation came from invariance properties of random walks and subsequently this question 
attracted the interested of many researchers and a large number of different approaches and generalizations were 
developed, each approach giving rise to a different stopping time which solves the same Skorokhod embedding 
problem — we do not even attempt to give a brief summary but refer to the surveys [27, 28, 21] which list many 
different approaches and discusses several applications. 

One of the earliest approaches after Skorokhod was initiated by Root in 1969 [30] who showed that if fj, is 
centered and has a second moment, then it is possible to find a a subset R of [0, oo] x [—00, 00], the so-called Root 
barrier, such that its first hitting by the time-space process (t,B t ), tr = inf {t > : (t,B t ) G R}, solves the 
Skorokhod embedding problem. In addition to this intuitive solution, Loynes [26] introduced extra conditions 
under which Root's barrier R is uniquely defined (as a subset of [0, 00] x [—00,00]) and Rost [34] proved that 
among all stopping times r which solve the Skorokhod embedding problem, Root's solution minimizes E [r 2 ] ; 
further he gave important generalizations and connections with potential theory, characterized the possible 
stopping distributions of Markov processes by introducing a filling scheme in continuous time and proved the 
existence of another barrier which leads to a hitting time which maximizes E [r 2 ] . Despite all these nice 
properties of Root's and Rost's approach they had until very recently the significant disadvantage that it was 
not known how to calculate the barriers — in fact, only for a handful of very simple target measures \x the 
explicit form of the barriers was known. 

One of the more recent motivations to study the Skorokhod embedding problem comes from mathematical 
finance where Skorokhod embeddings naturally appear if one is interested in model-independent prices of exotic 
derivatives, i.e. "extremal" solutions of the Skorokhod problem lead to robust lower and upper bounds for 
arbitrage-free prices of exotics — we refer to the survey's of Hobson [21] and Obloj [28]. Motivated by this, 
Dupire [14] gave in 2005 a presentation where he firstly pointed out that in the case of options on variance these 
extremal solutions are the ones of Root [30] resp. Rost [34] and, secondly, that these barriers can be calculated 
by solving nonlinear parabolic PDEs, namely solutions to free boundary problems. This was further developed 
by Carr-Lee [8] and in 2011 Cox- Wang [11] made this connection precise using the variational approach as 
developed in the 1970's by A. Bensoussan-J.L. Lions et. al, cf. [6] (and do much more — we refer the reader 
to [11]). This article is inspired by all these beautiful results and our goal is to study the connection of Root's 
barrier with the parabolic obstacle problem from the perspective of viscosity theory; more precisely, given two 
probability measures fi, v on (R, B (R)) and a function a : [0, 00) x R — > R we want to find (under appropriate 
conditions on cr, /i, v) the Root barrier R such that its first hitting time tr solves the following Skorokhod 
embedding problem 



(1) 



dX t = cr (t, Xt) dB t , Xq ~ fj,, 

X TR ~ v and X TR = (X tATB ) 4 is uniformly integrable. 



Key words and phrases. Skorokhod embedding problem, Root and Rost barrier, viscosity solutions of obstacle problems, reflected 
BSDEs, bounds on variance options. 
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(2) 



by identifying R with a continuous function u G C ([0, oo) x R,R) which is a viscosity solution of 

+ f M \.-x\v(dx) ~4 Au ) = on (0, oo) x R, 

u(0, .) = — | . — x| /x (dar) on R. 

and with the quadruple of process (X, Y, Z, K), solution to the reflected Forward-Backward SDE (for h (•) := 
- J* R I - - V\v (dy) , T > 0, (t, ar) G [0, T]xM and s G [i, T]) 

= u(0,X^. x ) +K t i x -Kl'*-£z t r '*dW r , 
Y*' x > h (Xl' x ) , t < s <T and jf (F/'* - h (X^ x )) d-ST*'* = 0. 
What we consider among the benefits of the viscosity approach are that it 

• covers the case when a is time dependent and only elliptic, i.e. no time-homogeneity or uniform ellipticity 
assumptions as in [11] are needed (especially the important case for model- independent bounds when 
a (t, x) — x is immediately in our setting without an extra transformation), 

• gives a direct link to reflected FBSDEs and their Feynman-Kac representation via classic work of El 
Karoui, Kapoudjian, Pardoux, Peng and Quenez[15], 

• uses continuous functions instead of functions with values in weighted Sobolev spaces and allows for 
short proofs of convergence of numerical schemes via the Barles-Souganidis method. 

Further, we introduce a generalization of Loynes' notion of regular barriers, so-called (/x, i/)-regular barriers 
which allows to establish a one-to-one correspondence between the problems (1) and (2). 

In Section 2 we recall several results from the literature and discuss the general Root embedding problem, 
in Section 3 we recall several results in viscosity theory and the connection with reflected FBSDES. Section 
4 contains our main result: the link between the embedding problem, PDE in viscosity sense and reflected 
FBSDE. In Section 5 apply the Barles-Souganidis [3] method in this context and we finish in Section 6 by 
recalling the relevance of Root's barrier for the derivation of bounds on variance options. 
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2. Convex order, potential functions, root barriers and root embedding 

Already for the case of Brownian motion one cannot hope to find Skorokhod embeddings for two arbitrary 
probability measures \x and u 1 i.e. Bq ~ /x, B T ~ v. Intuitively, a necessary condition is that the target 
distribution v should have a higher dispersion than the initial distribution \i. This can be expressed by a classic 
order relation on the set of probability measures, namely the so-called convex order. In the context of Root 
solution this order relation appears also naturally when one thinks in terms of the potential theoretic picture 
as developed in this context by Rost [31, 32, 33, 34], Chacon [10] and others. 

Definition 1. Let t £ N, m £ I and denote with Ad^ the set of probability measures on (R, B (R)) with a finite 
fc-th moment and mean m (i.e. /x £ -Mm iff L xfi (dx) = m and L \x\ k fj, (dx) < oo) and set Ai k = U m eR ^ti 
M = U eN u{o}- Mfe - 

Definition 2. Let m G R. A pair of probability measures (/x, v) G M 1 x A4 1 is in increasing convex 1 order, 
H < cx if 

(3) / / (x) \x (dx) < I f (x) v (dx) for every convex function / : R — > R 

Jr Jr 

(provided the integrals exist). We also say that two measures /z, v are in convex order if either /x < cx v or v < cx /i. 
Two real-valued random variables X, Y are in increasing convex order, denoted X < cx Y , if C (X) < cx C (Y). 

Remark 1. Assume /x G M.^, v G M.^ and [i < cx v. Then m = n, this follows from choosing / (x) = x and 
/ (x) = -x. 

The intuition being that /x < cx v signifies that a /i distributed random variable is less likely to take very large 
values than a v distributed random variable since extremal values of convex functions are obtained on intervals 
of the form (— oo, a) U (b, oo). As it turns out, an especially important case is given when (3) is applied with 
/ (.) = |. — m| resp. / (.) = — |. — m| for fixed m. 



^Sometimes this is also called the Choquet order. 
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Definition 3 (Potential function). Associate with every fi <E Ai 1 a continuous non-positive function u^ G 

C (R, (-00, 0]) given as 

Up(x) := - \x-y\fi(dy). 



We say that Up is the potential function of the probability measure fi. 

Example 1. Denote with S m the Dirac measure at m G R. Then u$ m (x) = — \x — m\ and it is easy to verify 
that Up < ug m for every fi G A^, resp. fi > cx o* m - 

Lemma 1. Let m G R and fi, v G A4„. TTien the following are equivalent 

(1) fi ^CI 

(2) - / (y - x) + fi (dy) > - f (y - x)+ v (dy),Vx G R, 

(3) u M (a;) > u v (x),\/x G R. 

Moreover, fi < cx v implies that f R xfi(dx) = f m xu(dx) and if fi,v G Ai^ then f ]gL x 2 fi(dx) < J" K x 2 i^ (<ia;). 

Proof. Follows from [35, Section 3. A]. □ 

Before we justify the name "potential function" let us note that the classic connection between positive 
[non-positive] Radon measures and second derivatives of concave [convex] functions reads in this setting as 

Proposition 1. There exists a one-to-one correspondence between potential functions and probability measures. 
More precisely: 

(1) Let fi G Ai^ for some m G R. Then Up (x) = — J R \x — y\ fi (dy) is a concave function satisfying 

Un (x) < u$ m (x) = — \x — ml and lim Up (x) — u$ m (x) = 

x— >±oo 

(2) Let u : R — > R 6e a concave function satisfying 

u (%) < —\% — m\ and lim u (x) + \x — m\ = 

/or some m G R. TTien £/iere exists a probability measure fi G A 7 ^ smc/i i/iai u is the potential function 
of pi, i.e. u = Up. 

Moreover, the connection is given by the relation fi ~ ^u" in the sense of distributions. 

Proof. See [20, Proposition 2.1] and [5, Proposition 4.1]. □ 

Remark 2. To explain the name "potential function" we briefly recall the classic potential theory for Markov 
processes (cf. [7, 27, 28] for details): consider a real- valued Brownian motion B and denote its semigroup 
operator with (P B ). The potential kernel is defined as U B = J °° P B dt, i.e. U B can be seen as a linear operator 
on the space of measures on (R, S(R)) by setting fiU B = fiP B dt which is of course nothing else than the 
occupation measure along Brownian trajectories started with Bq ~ fi. If fi is a signed measure with fi (R) = 
and finite first moment, then the Radon-Nikodym density with respect to the Lebesgue measure is given as 

Since (in dimension one) Brownian motion is recurrent, fiU B is infinite if fi is positive. However, the right hand 
side — f R \x — y\ fi (dy) is still well defined for every fi G Ai 1 and Chacon [10] demonstrated that this is indeed 
a very useful quantity in the study of Martingales and Markov processes with respect to hitting times. 

We summarize some properties of the potential functions of probability measures which we use throughout. 
The proofs are standard and can be found in [27, Proposition 2.3] and [28, Section 3.2]. 

Proposition 2. Let fi, v G Ai 1 . Then 

(1) fi < cx v iff fi > cv v iff Up (.) > u y (.), 

(2) Up is concave and Lipschitz continuous with Lipschitz constant equal to one, 

(3) if /j, (R) = m then Up (x) < u$ m (x) = — \x — m\, 

(4) if v G M 1 then V6 G R we have 

Up\[b,oo) = Wi>|[b,oo) iff M|(6,oo)= H(h,oo), 

(5) if pi n G Ai 1 , fi n =^n fi- weakly 2 and Up n converges for some point xq G R then a Mn converges for all 
x G R and 3c > s.th lim„_ ! . 00 Up n = u u + c. Further, if 3 v G Ai 1 such that u Un > u v then fi G Ai 1 
and c = 0. Conversely, if u Url (x) —> n (x) V.t G R for some fi„, fi G Ai 1 then fi n fi weakly, 

(6) if fi < cx v then 1 &m\ x \-+ oa (u^ (x) - u v (x)) = 0. 

(7) u a is almost everywhere differentiable on R and its right derivative at x equals 1 — 2fi ((—00, x]) and its 
left derivative equals 1 — 2fi ((— 00, x)). 



2 The symbol "=>" denotes weak convergence of probability measures. 
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Since we want to embed a measure via a Markov process, the potential function which is given by the image 
measure of this Markov process is a fundamental quantity. 

Definition 4. If X is a real-valued, continuous Markov process carried on a probability space (f2, F, (Ft) , P) 
with P o XfT 1 £ M 1 for every t > then we call 

u x (t,x) := u PoX -i(x) 

the potential function of X (at time t). 

As shown by Chacon [10, Lemma 3.1a], potential functions give an intuitive characterization of Martingales. 

Proposition 3. Let X be a real-valued, adapted process. Then X is a martingale iff for every pair T\,Ti of 
bounded stopping times with T\ < T2, 

ux (n, ■) > ux (t 2 , ■) ■ 

2.1. Root's solution and classic results on barriers. Our aim is to provide explicit solutions of the Sko- 
rokhod embedding problem not only for Brownian motion but for diffusion martingales. 

Assumption 1. (f2, F, (Ft) , P) is a filtered probability space that satisfies the usual conditions and carries a 
standard Brownian motion B. Assume that a £ C ([0, oo) x K,R) is Lipschitz in space and of linear growth, 
both uniformly in time, that is 

. . [<7 (t, x) — a (t, y)\ \a(t,x)\ 

(4) o~Lip '■= sup sup j j < oo and g^q := sup sup — : — - < oo. 

t£[0,oo) x^y \ x — U\ te[0,oo)zGR 1 + \ x \ 

Further, let ra 6 1, (p,, v) 6 M^, p < C x v, ond Xq ~ p, Xq independent of the a-algebra generated by the 
Brownian motion B and assume that with probability one a 2 (t,X t ) > Vt > where X = (X t ) denotes the 
unique strong solution of the SDE (6), i.e. the real-valued, local martingale X such that 



Further, assume that for every t > 0, X t has a continuous (but not necessarilty smooth) density with respect to 
Lebesgue measure. 

Definition 5 (Skorokhod embedding problem). Let (a,fj,,v) and (f2, J 7 , (.Ft) , (X t ) ,P) be as in Assumption 1. 
If there exists a stopping time r (wrt the filtration (Ft)) such that 

X T ~ v, X T = (Af AT ) t>0 is uniformly integrable (u.i.) 

then we say that r solves the Skorokhod embedding problem given by (a,/j,,v). We denote with S{a,fjL,v) the 
(possibly empty) set of stopping times r which solve the Skorokhod embedding problem given by (a, /i, v). 

Remark 3. Unless stated otherwise, all stopping times are stopping times in the filtration (Ft). 

As pointed out in the introduction there are many different solutions of the embedding problem, that is in 
general the set S (a, (i, v) can contain much more than one stopping time (in the survey article of Obloj [28] at 
least 20 different solutions are given for embedding laws in Brownian motion, i.e. (a, fi, v) = (1, Sq, v)). In terms 
of applications a natural requirement is to consider stopping times which are minimal. 

Definition 6. A stopping time r is called minimal for the process X if for every other stopping time p such 
that p < t and X T ~ X p it follows that p = r. 

This article is focused on the approach initiated by Root in 1968 [30] and further developed by Rost, Loynes, 
Kiefer et. al [34, 32, 31, 25, 26] of solving Skorokhod embeddings by stopping times which are the first entrance 
times of the time-space process (t, X t ) t>Q into a closed subset R of [0, oo] x [— oo, oo]. 

Definition 7. A closed subset R of [0, oo] x [— oo, oo] is a Root barrier R if 

(1) x) e R implies (t + r, x) E R Vr > 0, 

(2) (+oo, x) £ R V.t e [-00,00], 

(3) (i,±oo) e RVte [0,+oo]. 

We denote by 1Z the set of all Root barriers R. Given R £ 1Z, its barrier function fp : [—00, 00] — > [0, 00] is 



Barrier functions have several nice properties such as being lower semi-continuous and that (fn (x) ,x)<ER 
for any i£l (see [26, Proposition 3] for more properties). 



(5) 




^0 



defined as 



/fl (x) := inf {t > : (t, x) £ R} , x £ [-00, 00] . 
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Definition 8. We denote with TZ (er, fi, v) as the (possibly empty) subset of TZ of all R G TZ (er, fi, v) such that 
tr <E S (er, /i, v) where tr is the first hitting time of R by (t, X t ), i.e. 

T R = mi{t>0:(t,X t )eR}. 

We also say that R G 1Z (er, /i, v) embeds the law /z into v (via X) and call tr the Root stopping time. 

Before we discuss conditions on er, /i, v which imply the existence of a Root barrier, i.e. that 1Z (er, fi, v) ^ 0, 
we note that simple examples show that the Root barrier is in general not unique, i.e. \1Z (er, fi, v) > 1 already 
for generic situations. 

Example 2. The set R = [0,oo] x [1, oo] 1J [0, oo] x [— oo, — 1] is an element of 1Z (l, So, | (S-i + Si)), i.e. R 
embeds the sum of two Diracs into a standard Brownian motion but obviously one can define R differently, 
in fact any other Root barrier with a barrier function which coincides with Jr on [—1, 1] is also an element of 
TZ (1,60,1(5.! + <*!)). 

This non-uniqueness problem was in the Brownian case resolved by Loynes [26] in 1970 who introduced the 
notion of regular Root barriers. 

Definition 9 (Loynes [20] regularity). Let R G TZ. We say that R resp. its barrier function Jr is Loynes-regular 
if Jr vanishes outside the interval [x R , x^\ , where x+ and x R are the first positive resp. first negative zeros of 

Jr. Given Q, R G TZ we say that Q, R are Loynes- equivalent if Jq = Jr on x_, x+ and \x R , x R 



-J 



Remark 4. If Q,R are Loynes-equivalent then x+ = x+ and 



x_ 



We advise the reader to be suspicious about the word "regular" since Loynes-regular Root barriers can have 
spikes (see examples in Section 2.2). For the case (er, n, u) = (1, So, v) existence and uniqueness of Loynes-regular 
Root barriers are standard. 

Theorem 1 ([30, 26]). For every v G A4q the set TZ (1, So, v) is not empty and there exists exactly one element in 
TZ (1, So, v) which is Loynes-regular. The second moment of the Root stopping time, E [t^] for R G TZ (1, <5o, v), 
is minimal among all elements of S (1, So, v). 

An advantage of Root's solution besides its intuitive appeal (the stopping time being a hitting time) the 
property that tr minimizes E [r 2 ] among all solutions of the Skorokhod embedding problem, i.e. in our notation 
E [tJ,] = inf E [r 2 ] where the inf is taken over all r G S (a, fi, v). This was proven in Rost [34] (as well as in 
[11]) and is the very reason why Root's solution is important for financial applications, see Section 6. 

2.2. Regularity beyond Loynes. We now extend the discussion on Root barriers to our more general setting 
of Root embeddings for Ito-martingales with random initial distribution, i.e. we discuss TZ (a, fx, v). To that end 
we need to adapt several results of Loynes [26]. First, note that already in the Brownian case (a = 1) we need a 
modification of the notion of Loynes-regularity, a notion that is tailor made for the Dirac as initial distribution 
as the example below shows. 

Example 3. Let \x = \ {Si + S-2) and v = j (S3 + 61 + <5_i + £-3). It is easy to see by symmetry properties 
of the Brownian motion that for a = b = the barrier 

Q atb = [0, 00] x [3, 00] (J [a, 00} x {1} (J [b, 00} x {-1} (J [0, 00] x [-00, -3] |J {+00} x [-00, +00] 

is an element of TZ (1, fj,, v), as is 

R = [0, 00] x [3, 00] (J [0, 00] x [1, -1] (J [0, 00] x [-00, -3] |J {+00} x [-00, +00] . 

However, neither is Loynes-regular and in fact there cannot exist a Loynes-regular barrier in TZ (1, /1, v). (Assume 
Qa,b G 7Z(l,fj,, v), then a, b > otherwise it would not be Loynes regular; now note that Qo,o G TZ(1, fi,is), 
hence every other Q a ,b puts under v more mass on the point 3 than the required \ since the geometry of Q a ,b 
implies that only more trajectories can hit the line [0, 00] x {3} than in the case a = b = 0; further, every 
element of TZ (1, fi, v) must coincide with Qo,o on [0, 00] x [1, 3] U [0, 00] x [—3, —1]). 

Motivated by the above we introduce the notion of (fx, v)-regular barriers. 

Definition 10. For /j, < cx v define 

N^ v := {x G (-00, +00) : Ufi (x) = u u (x)} U {±00} and M^ v = [0, +00] x N^ v . 

A Root barrier R is (fx,v) -regular if R = RUj\f fl ' L ' [or equivalently if fn.(x) = Vx G N^- v ]. Denote with 
TZreg (cj Mi v ) the subset of TZ (a, [i, v) of (/i, z/)-regular Root barriers. Further, two Root barriers R, Q are (jx, v)- 
equivalent if 3 R \ ([0, 00] x (N^ u )°) =Q\ ([0, 00] x (N^)°) [or equivalently if f R (x) = f Q (x) Vx G (N^) c ]. 



^We denote with A the closure and with A" the interior of a given set A. 
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If A M and F v are the cumulative distribution functions of [i and v respectively, then the set N^'" coincides 
with the set {x : (x) = F v (x)} U {±00}. The next Lemma shows that in the case a = 1 and \i = So the 
(/x, j/)-regularity of a Root barrier is equivalent to Loynes regularity. 

Lemma 2. Let v G M\ and R G 1Z (1, <5o ; v). Then R is Loynes-regular iff R is (So ,v) -regular. 

Proof. Assume R to be Loynes-regular, then fu (x) = for Vx ^ (x^,x^J where x^,x^ are as in Definition 
9. Since the embedding is possible, i.e. X TR ~ v, then via the regularity of the barrier R one has that 
^([x^, #+]) = 1, i.e. v puts no mass outside [x^,x+] (by definition so does So). In other words outside 
[x^, x+] both measures v and i5o are the same which by point (4) of Proposition 2 implies that ug (x) = u v (x) 
for x ^ [x^,x+] and by continuity of the potential functions for any x ^ (x;?,x+). We conclude that R is 
(So, ^)-regular. 

For the inverse direction, just remark that being R a (So, ;/)-regular barrier we have due to So < CI v that 
N s °' v n K has the form TV 50 -" n R = R\ (a, b) for some a < < 6, in other words f R (x) = for any x £ (a, b). 
Notice that the convex order relation (with So) implies that if u v (c) — for some c > then u v (x) = for any 
x > c (and corresponding condition for the case c < 0), this implies that a and b are the first negative resp. 
positive zero of fu . Hence R is Loynes-regular. □ 

Definition 10 allows to extend Loynes' [26, Theorem 1] to the present setting, i.e. if there exists an embedding 
with a Root barrier then there also exists an embedding with a (/i, j/)-regular Root barrier. 

Theorem 2. Let /1, v G M 2 , [i < cx v and let X be a real-valued continuous square integrable local Martingale 
X with a quadratic variation that is strictly increasing and Xo ~ (J.. Assume additionally that there exists 
Q G 1Z s.th. X TQ ~ v and X TQ is u.L. Then Q is (n,v) -equivalent to a unique (fj,, v) -regular barrier R G 1Z 
s.th. X TR ~ v and X TR is u.L; moreover E [[X] TR ] < 00. 

We prepare the proof with a lemma concerning the union and intersection of Root barriers. 

Lemma 3. Take v G Ai 1 and let Y be a process with continuous paths. Assume there exists Q,R G 1Z with 
first hitting times tq and tr respectively and assume that Y TQ ~ v and Y TR ~ v . Then Q U R,Q D R G 1Z and 
both generate the same law v with respective stopping times min {tq,tr} and max {tq,tr}. 

Proof. This proof is a straightforward modification of [20, Proposition 4 ], we fill in some details and sketch 
the proof for the case of Q U R: let fu and /q be the barrier functions of R and Q and define the set 
K := {x : fq (x) < fu (x)}. Since both barriers (with their respective stopping times) generate v we have that 
P [Y TQ G K] = P [Y TR G K] iff 

P [Y TQ G K, Y Ta eK]+¥ [Y TQ G K, Y TR G K c ] = P [Y TQ G K, Y TR G K] + P [Y TQ G K c , Y Tr G K] 

iff P [Y TQ G X, F TJ7 G K c ] = P [y TQ G K c , Y TR G A'] . A quick analysis shows that for a trajectory Y (w) s.t. Y TR (uj) 
ends in K it is impossible that Y TQ (w) G A c , in other words, since the paths of Y are continuous and in the set 
[0, 00) x K the barrier Q is to the left of R then any path that hits Q must hit R before (hitting Q) and hence 
P [Y TQ G A c , Y TR G K] = = P [y TQ G A, F TJJ G K c ] . 

Finally, by decomposing the probability space in disjoint events (based on A and A c ) and using the just 
obtained result one concludes P [{Y TQ G A, Y TR G A} U {Y TQ G A c , Y~ TH G A c }] = 1. The rest of the proof 
follows as in [26]. □ 

We can now give the proof of Theorem 2. 

Proof of Theorem 2. To see that Q is (/i, ^-equivalent to a (n, z^-regular barrier note that since u v < (and 
we embed by assumption) the continuous time-space process (t A tq, X t /\ TQ ) does not enter [0, 00] x N^' v , hence 
R := QU M^' v is a also an element of 1Z (c, /1, v), and moreover an element of 7Z reg (cr, /i, v). 

To see the uniqueness, i.e. \7Z reg (a, v)\ = 1 suppose that there are two (/1, z^)-regular barriers B,C, each 
embedding ^ (via A) with u.i. stopping times tb and tq respectively. Then r = B U C also embeds v with 
stopping time 7 = min {tb,tc} (see Lemma 3); 7 is also u.i. Furthermore, since A is a local martingale, 

E [[A] TB ] = E [A T 2 J = E [A 7 2 ] = E [[A] 7 ] . 

Since the quadratic variation t i-> [A] t is strictly increasing and 7 < tb this already implies [A] 7 = [A] TB . 
The strict monotonicity of 1 1-> [X] t then implies 7 = tb- The same arguments applies to rc, hence we conclude 
that = 7 = rc 1 . 

To see that B,C and BU C are the same (outside A/" M:iy since in TV^'" this must hold) we argue as in the 
proof of Lemma 2 in [26] by showing that if B ^ T then also tb ^ r 7 . Note that B / T implies that the 
existence of xo £l such that xo ^ iV 4 '" and (wlog) (xo) < fr (^0) where fr (xo) 7^ 0. Since the barriers are 
regular and lower semi-continuous functions attain their minimum on a compact set, 3e > s.th fr (x) > e on 
some neighborhood of xo- Moreover e can be chosen in such a way that 2e < fr (xo) — /b (#0)1 implying that a 
<5 > exists s.th fr (x) > /r (xo) — e > /b (xo) + e for any x in the set {x : |x — Xo| < S}. Since the embedding 
via X is possible, there exists a positive probability that the time-space process (i, X t ) t>o hits the line segment 
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[fs (xo) j Ib (xq) + e] x {xq} while hitting T only after time fg (xo) + £■ This implies that P [tb 7^ Tr] > and 
the result follows. Since Q 6 72. (a, [1, u), E [[X] TQ \ < 00 and hence by construction also E [[X] TR ] < 00. □ 

2.3. The general Root embedding and optimality. We give a generalization of Theorem 1 showing that 
1Z (a, /1, v) is not empty, i.e. embedding v E M. 2 via the solution process of the SDE 

(6) dX t =a(t,X t )dB t , X 

Root's original construction (for Brownian motion, i.e. a = 1 and [i = So) relies only on the continuity of the 
Brownian trajectories, the strong Markov property, that trajectories exit compact intervals in finite time and 
the fact that the process has no interval of constancy - hence (as remarked several times in the literature) the 
same construction generalizes. Here we give only the results, complete proofs are postponed to Appendix 7. It 
turns out to be natural for the proof to distinguish between the cases [X] ^ = 00 and when X is a geometric 
Brownian motion (which implies only [In A] oo = 00). 

2.3.1. Existence and uniqueness of a Root barrier if [X] x = 00. We now state a general result on the existence 
of the Root embedding via the solution of (6) (the proof is given in the appendix). 

Theorem 3 (Root Embedding). Let (a,[i,i>) be such that Assumption 1 is fulfilled and that [X] ^ = 00 a.s. 
Then \lZ reg (cr, fi, v)\ = 1. Moreover, 

E [[X] TR ] = J x z v(dx) < 00 and E [tr] < 00 Vi? G TZ(a,fi, v) . 

Remark 5. One can also apply the results of Rost and Cox- Wang [11, 34] to show that tr minimizes E [/ (r)] 
among all r G S (a, /1, v) where / can be a convex, increasing function with / (0) = which further fulfills some 
properties which depend in general on a. 

2.3.2. Existence and Uniqueness for geometric Brownian motion (a (t, x) = x). The case when X is a Geometric 
Brownian motion (GBM henceforth), i.e. a (t, x) = x, is not covered by Theorem 3 since although E [[X] = 00 
we have [X] ^ < 00 a.s. (otherwise one could write X as a time-changed Brownian motion via Dambis/Dubins- 
Schwarz). However, this case is especially important for the applications in finance (see Section 6). 

Assumption 2. (//, v) € M 2 , [i < cx v and there exists an e > such that 

supp (fx) C (0, +00) and supp (y) C [e, +00) . 

The assumption that e is strictly greater than is needed for our arguments (see Remark 17 in the appendix) 
and does not cover the case v ([0, 00)) = l,zy({0}) = 0. (However, note that one cannot hope for Root 
embeddings where the target measure has an atom at 0). We now state the embedding result for this situation 
(again, the proof is given in the appendix). 

Theorem 4 (Root Embedding via Geometric Brownian motion). Let a (x) = x and Assumption 2 hold. Then 
\R, reg (o-,\i,v)\ — 1. Moreover 

E [[ X \ r R ] = / x 2 fi (dx) < 00, [InX] TR = t r VREK (cr, fi, v) 

Further, for every other stopping time t € S (cr, /1, v), 

E[/(r fi )]<E[/(r)] 

for every function f : [0, 00) — > [0, 00) which is convex, f (0) — and f has a bounded right-derivative. 

3. Background on viscosity solutions and the obstacle problem 

We briefly recall the definition and some basic stability properties of viscosity solutions in the parabolic 
setting (for a detailed exposition we refer to the User's guide [12]) with some additional lemmas which will 
become useful in the subsequent sections. We then discuss the connection to reflected FBSDEs. 

3.1. Sub- and supersolutions, properness and semirelaxed limits. 

Definition 11. Let O be a locally compact subset of K and denote Ot = (0, T) x O for T G (0, 00]. Consider 
a function u : Ot — > K and define for (s, z) G Ot the parabolic superjet Vq + u (s, z) as the set of triples 
(a,p, m) G K x R x R which fulfill 

u (t, x) < u (s, z) + a(t — s) + (p, x — z) 

(7) +i (m (x - z) , x - z) + o (\t -s\ + \x- z\ 2 ^j as O t 3 (*, x) -» (s, z) 
Similarly we define the parabolic subjet Vq~u (s, z) such that Vq~u = —P 2 '+ (—it). 
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Example 4. If a function u £ C 12 (Ot,R) then 



du 


du 






du 


du 








(t,x) = 



du 



\(t,x) 



V 2 ' + u(t,x) U ^\(t,x),^\(t,x),m\ \( t , x y.m> J^L| (t 

hence 

' du du du 
dt' dx' d 2 x t 

(If u is only differentiable from the left or right, a similar statement holds and is important for the proof of our 
main result, see Lemma 5). 

Remark 6. An equivalent definition is to consider the extrema of u — ip when tp is taken from an appropriate 
set of smooth test functions, see [12]. 

Definition 12. A function F:G> T xRxRxRxR^Ris proper if V (t, x, a, p) £ Ot xlxl 

F (t, x, r, a,p, m) < F (t, x, s, a,p, m!) Vm > m! , s > r. 

Denote the real-valued, upper semicontinuous functions on Ot with USC (Ot) and the lower semicontinuous 
functions with LSC (Ot)- A subsolution of the (forward problem) 

F (t,x,u,d t u, Du, D 2 u) = 



(8) 

is a function u £ USC (Ot) such that 

O 



u(0, .) = u (.) 



F (t,x,u,a,p,rn) < for (i, x) £ Ot and (a,p, m) £ Pq + u (i, x) 



w (0, .) < u (■) on 

The definition of a supersolution follows by replacing USC (Ot) by LSC (Ot), V@ + by V@~ and < by > . If u 
is a supersolution of (8) then we also say that F (t, x, u, d t u,Du,D 2 u) > 0, u(0, 

x ) ^ uq (x) holds in viscosity 

sense (similar for subsolutions). Similarly we call a function v a viscosity (sub-, super-) solution of the backward 
problem 

G(t,x,v,d t v,Dv,D 2 v) = 

v(T,.) = v T (.) 



(9) 
if 



G (t,x,v,a,p,m) < for (f, i) e Or and (a,p, m) £ Pq + v (t, x) 
v (0, .) < v T (■) on O. 

The next lemma shows that finding a viscosity solution to (9) is equivalent to finding one to (8). 

Lemma 4. Let T < oo. Then u is a viscosity solution of the forward problem (8) iff v (t, x) := u(T — t, x) is 
a viscosity solution of the backward problem (9) with Vt (•) = Uo (.) and 

G (t, x, v, a,p, m) = F (T — t,x, v, —a,p, m) . 

Proof. Assume u is a viscosity solution. Let ip £ C°° (Or,R) and assume (t,x) i-> v(t,x) — ip(t,x) attains 
a local maximum at (t,x) £ Ot- Then (T — t,x) is a local maximum of (t,x) i— > u(t,x) — ip(t,x) with 
ip (t, x) = <p (T — t,x), hence 

F (; ; U, d t p>, Dip, D 2 p) | (T _ U) = G (; -,v, -8t<p, Dip, D 2 p) | (t>) 

which is by Remark 6 sufficient to see that v is a viscosity solution of (9). The other implication follows 
similarly. □ 

A strong feature of viscosity solutions, which is also the key for the proof of our main theorem, is that 
they are very robust under perturbations. The type of reasoning presented below is classic and called the 
Barles-Perthame method of semi-relaxed limits. 

Definition 13. Let A C K m and (g n ) n a sequence of functions g n : A — > R which is locally uniformly bounded. 
We define g, g : A — > R as 

g(p) = limsup g n (q) . 

q—tp,n— too 

g(p) = liminf g n (q) . 
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Proposition 4. Let {u n ) n C U SC {Ot), O a locally compact subset o/R, {F n ) n a sequence of functions 

F n : Ot xlxlxl^l 

such that each u n is a subsolution of 

F n (t,x,v,d t v,D 2 v) = 0. 
Further, assume (u n ) and (F n ) are locally uniformly bounded. Then u is a subsolution of 

F (t, x, v, d t v, D 2 v) = on O. 

The analogous statement holds for a sequence of LSC {Ot) functions which are supersolutions. Further, if 
u = u then the convergence of (u n ) to u = u is locally uniform. 

Proof. This is simply a restatement for parabolic PDEs of Proposition 4.3, Lemma 6.1 found in the User's guide 
[12] combined with their subsequent remarks: the definition of u as limsup,^^ ^ x ) u n (s, y) guarantees 
the existence of sequences 

rij — > oo, (tj,Xj) G Ot s.t. (tj,Xj) — > (t,x) and u nj (tj,Xj) — > u 

and it is clear that for every sequence Ot 9 {tj, x j) converging to (t,x) one has lim sup j u„ ^ (tj,Xj) < u. This 
allows to apply [12, proposition 4.3 ] which implies that V(a,p, X) G T j2 ' + u (t, x) there exists a sequence 

(tj,Xj) C Or, (pj,m,j) G Vo + u nj {tj,x 3 ) s.t. (tj,Xj,u nj (xj) ,aj,rrij) (t, x,u(x) , a,m) G V^ + u{t, x) . 

However, each u nj is a subsolution so that for every j it holds that 

i* •".;•/•';• m j) < o 
and hence from the definition of F_ it follows that 

F_ (t, x, a,p, m) < 0. 

The proof for subsolutions follows the same argument. Finally, if u = u =: u then u must be continuous and 
the convergence u n — > u must be uniform (as seen by [12, Remark 6.4]). □ 

3.2. Reflected FBSDEs and viscosity solutions of obstacle problems. Our goal is to show that knowing 
the Root barrier is equivalent to knowing the viscosity solution of a certain obstacle PDE, that is we are interested 
in the case F (t, x, r, a,p, to) = min (^r — h (x) , a — a y^- roj where the initial condition uq and the barrier h are 

given as the potential function of the initial and target measures and a is continuous (we discuss the motivation 
for looking for such a PDE solution at the beginning of Section 4). 

Definition 14. Let T € (0, oo], u : K R and h, a G [0, T) x R -> R. Denote with V T (cr, u , h) the (possibly 
empty) set of viscosity solutions 



(10) 



(u(t,x)-h(x),(d t -^-A)u(t,x)\ = 0, (t, x) G (0, T) 

u(0, .) = u (x),xeR. 



which are of linear growth uniformly in t, i.e. V (t, x) G [0, T) x R, sup^ t x ^ e \^ t)xR < 00 ■ ^ e re ^ ev t° this 

PDE as the obstacle problem on (0,T) x R with barrier h. 

Remark 7. F (t, x, r, a,p, to) = min (r — h [x) , a — a Y^- rnj is proper and moreover the barrier h introduces 
an asymmetry in the definition of sub- and supersolutions: in this case Definition 12 says that u is a subsolution 
of (10) if u (0, .) < uo (.) and for every (t, x) G (0, T) x R such that u (t, x) > h (x) we have 

a 2 (t x) 

(11) a ^m<0V(a,p,m)ePg' + ti(t,x). 

Whereas u is a supersolution if u (0, .) > uq (.) and V (t, x) G (0, T) x R 



cr 2 (t, x) 



(12) a — — -m > V (a,p, to) G V^' u(t,x) and u (t, x) > h (x) . 

We now discuss the abstract problem of existence and uniqueness of a viscosity solution to (2). First we recall 
the definition of a subclass of reflected forward-backward stochastic equation (RFBSDE) ([ 1 5, Section 8] applied 
with / = 0). Therefore fix T G (0,oo) and for each t G [0, T] denote with {F*,t < s < T} the natural filtration 
of the Brownian motion {W s — Wt, t < s < T}, augmented by the P-null sets of T . Take (t, x) G [0, T] x R and 
define for s G [t, T] the forward-backward dynamics whose solution is the ^-adapted quadruple (X, Y, Z, K) 

Xt' x = x + J t s a(T-r, X l r x ) dW r , 

(13) { = //„(A7'') • Ay /C j'J /T,"'/U ,. 
Y^ x > h (X*> x ) , t < s < T and jf (Y*> x - h (X*> s )) dK^ = 0, 
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where (K l s x ) sG r t is an increasing and continuous process verifying K t ,x = 0. In the above equation the 
meaning of the processes X and Y is clear and a rough interpretation of the processes Z and K is a follows: Z 
guides the evolution of Y via the Ito integral so that Y can hit the random variable uo (Xt) at horizon time 
T. Note that both processes Y, Z are (J-t)-adapted and so they do not "see the future", nonetheless, uo (Xt) is 
attained at time t = T. The process K ensures that Y does not go below the barrier h: K pushes Y upwards 
whenever Y touches and tries to go below the barrier h, else it remains inactive (that is constant). 

Proposition 5 (Stochastic representation for the obstacle PDE). Let uq, h G C (R,R),cr G C ([0, oo) x R,R) ; 
h (•) < uq (.) and assume that for every t G [0, oo), a (t, .) , h (.) , uq (.) are Lipschitz continuous (a uniformly in 
t). Moreover, assume that 

K (a?) | + \<r(t,x)\ + \h(x)\ 



sup sup 

te[o,oo) iei 



< oo. 



l + \x\ 

Then there exists a G C7([0,oo) x R, R) which is a viscosity solution of 

lin (v (t, x)-h (x) , (d t - £ A) v (t, x)) 

v(0,x) 

Moreover, for every T G (0, oo) 

(1) there exists V (t,x) G [0, T] xR a unique quadruple (Xl' x , Y*> x , Z\' x , K\' x ) s <=u,t] of {J 7 ^} -progressively 



(14) 



0, (t,x) G (0,oo) x 
wo (a;) , x G R. 



measurable processes which satisfies for some ct > 
sup {\Xl> x \ 2 + \Y s t ' x \ 2 } + 



sup E 

te[o,T] 



»e[t,T] 



zi 



l dr + Kl 



< C T (1 + |x|) 



a«d this quadruple fulfills (13). Moreover, V(s,x) G [i, T] x R i/ie mapping (s,x) i-> (-Xj'" 1 , y s ) is 



-a.s. jointly continuous and in particular (xl' x , Y^'^, Z t t,,J: , _fC. 

(2) (t,x) i — ^ Y^,Zl' X is the unique element of Vt (o~,uq, h) 

(3) Woo|[0,T]xe(^,x) = 1^ 



(i.e. deterministic), 



T-t,: 

T-t 



min ( v (t, x) — h (x) , ( — dt — %- A ) u (i, x) 

w (T, x) 



Proof. Existence & uniqueness in [0, T] x R, T < oo; Time reversion of (2) (see (9) and Lemma 4) shows that 
it is enough to deal with 

0, (t,x) G Ot = (0,T) x R 
Mo (x) , X G R. 

where a(t,x) = o~(T — t, x). Continuity, existence and uniqueness of the viscosity solution v follows from 
Lemma 8.4, Theorems 8.5 and 8.6 in [15] respectively. The linear growth of u in its spatial variable follows from 
standard manipulations for RBSDEs. [15, Proposition 3.5] applied to the RFBSDE setting above (i.e. (13) with 
a replaced by a) yields the existence of a constant Ut > such that V (t, x) G [0, T] x R 



Y, 



t,x\ 2 



< E 



sup | Y; 

se[t,T] 



< k T (E[|uo (Xt X )\ 2 ] +E 



sup \h + (s,Xl 

s£[t,T] 



<c T (l 



with h + := max{0,/i} and where the last inequality follows from the linear growth assumptions on uq and h 



E 



I vt,x | 2 
T 



] < or (1 



(see e.g. equation (4.6) of [16]). The 



along with standard SDE estimates: sup te [ T \ 
solution to (2) follows via our Lemma 4. 

Existence & uniqueness in [0, oo) x R: Above gives a unique solution for every finite T > 0. First note 



that for T,T' > 0, Y^Z t and Y^Zt coincide on [0, T A T') x R. By the comparison result in the appendix, 



rT'-t 



Theorem 9, we can define (t, x) as Y^_f ,x where T can be arbitrary chosen subject to T > t. Since the 
definition of a viscosity solution just relies on local properties, is a viscosity solution of (14). □ 

Remark 8. In the above link between the solutions of (13) and (the time reversed version) of (10), (14) one 
can in general not expect that the solutions of (14) exist in a classical sense. The following formal argument 
gives at least an intuition why RFBSDE and obstacle PDEs are in a similar relation as SDEs and linear PDEs: 
suppose a sufficiently regular solution v of (14) exists. Via Ito's formula it follows that: 

—2 



Y^ = v (s,Xi< x ) , ZY = (vS7 x v) (s,X^) and K* a 
The last condition in (13) then reads as 

J (Y^' x -h(X t r ' x ))dK t r ' x = 0iS J (v-h)[ 
and the rhs explains the form of the PDE fulfilled by v. 



-dtV — —Av 



,Xt' x )dr. 



-Av 



(r, Xl' x ) dr = Q 
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Remark 9. Proposition 5 applies as well to moving barriers (i.e. h could also depend on time) however we only 
need barriers constant in time. Note also that the above shows {i>oo|[o,t]xr} = Vt (c, u^jU^) but it does not 
show that {«oo} = "^oo {o',u fll u l/ ) 1 viz. we only claim that Voc is a viscosity solution of the PDE (14) on the 
unbounded time domain but not that Vac is of linear growth uniformly in time (however, we this will be a simple 
consequence of Theorem 5) . 

The existence proof of (13) via penalization in [15] gives a lower bound for the first time the solution touches 
the barrier. This is intuitive clear but since it is useful for numerics we give the full proof. 

Corollary 1. Let Woo be the solution of (14) as given by Theorem 5 and let w to be the unique viscosity 
solution 4 " of 

(15) f (d t - ^ ^ w° (t,x) = 0, on (0,T)xR 

\ w°(0,-) = «„(•) 

Then 

w° (t, x) < Woo (t, x) V (t, x) G [0, T] x R. 

In other words, Voo hits the barrier h at time t only after w° has hit it at some time < s < t. Furthermore, 
if a is uniformly parabolic and bounded then w° G C 1,2 ((0,T) x R,R) and can be written as a fundamental 
solution. 

Proof. The statement follows by noting that Woo is a viscosity solution of (14), hence a supersolution of the 
heat equation (15). Applying a standard comparison result immediately shows w° < Wqo. Another way to see 
this without resorting to comparison results and which gives another construction of a solutions to the obstacle 
problem is the following: the existence proof for (13) was done in [15] via penalization arguments. For neH 
consider the penalized obstacle problem 

d t w n (t, x) - iff 2 (t, x) Aw n (t, x)=n {w n (t, x) - h {x))~ , (t, x)eO T , 

w n (0, x) = u (x) , iei 

where a~ := — min {a, 0} for any a G R implying (w n (t, x) — h (x)) := — min {w n (t, x) — ft. (x) , 0} . Intuitively, 
the solution of the penalized equation approximates that of (2) as n — > oo since the term n (w n — h) ~ blows 
up for the points which do not satisfy the barrier condition and so in order for a well behaved solution to exist 
in the limit the barrier condition has to be verified. To shorten the discussion, along the proof it is shown the 
sequence (w n ) „ £ n converges and that w n j" u (as n — > oo), i.e. the sequence is monotonic increasing (see [15, 
Theorem 8.5]). It follows that w° (t, x) < u (t, x) for all (t, x) G [0, T] xK with w° the solution to the penalized 
PDE for n = 0, equation (15). 

Standard results (see for example [18, Chapter 1, Section 7, Theorem 12]) guarantee that under the stronger 
assumption of boundedness and uniform parabolicity of a there exists a smooth solution and since every smooth 
solution is a viscosity solution the statement follows. □ 



4. Root's barrier and the parabolic obstacle problem 

We are now in position to prove the main results: a one-to-one correspondence between regular Root barriers 
and viscosity solutions. 

4.1. From the Root barrier to the obstacle problem. As pointed out in Section 2, potential functions 
are a powerful way to keep track of the evolution of the distribution of a local martingale. Let R G 1Z (a, fi, v) 
and consider the potential function of the stopped process X TR , namely u (t, x) = ux^r (t, x) = E [\Xl R — x|]. 
Tanaka's formula immediately gives 

u (t, x) = u M (x) - E [L x tfXTR ] 

where (Lf ) t denotes the local time of (X t ) t , i.e. we see for every fixed x G R that 1 1->- u (t, x) decreases until 
we have substracted enough such that it matches u v (x). While one cannot hope that u is a classical PDE 
solution (due to the kinks that appear when u touches the barrier u v ) we give a formal derivation for the PDE 
satisfied by u which then motivates our proof of Theorem 5. For brevity also assume a = 1 (i.e. we embed 
into Brownian motion). Consider a sequence (/")„ C C°°(R,R) such that f n (.) — > n ^oo \-\ uniformly, hence 
V/" (.) sgn (.) and A/" (.) ^ n So (.)• Set 

(16) u n {t,x):=-nf n {XV'-x)]. 



This equation is a standard parabolic PDE whose unique solution is known to exist under the assumptions of Theorem 5, see 
e.g. the proof of Theorem 8.5 in [15] or think of the RBSDE (13) when h = — oo (implying that Kt = 0). 
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and since /" is smooth, Ait™ (t, x) = — E [A/™ (X^ R — x)] and from Ito's formula applied to the right hand side 
of (16) it follows that 

(17) u n (t, x) = - [ f n (y-x)fi (dy) - i f E [A/" (X r - x) l r < TR ] dr. 

Jr 1 Jo 

Hence dtu n (t, x) = — ^E [A/™ (X r — x) lt< Tr ] an d putting the above together shows 
d t u n - ^Au 11 = ^E [A/' 1 (X T t R -x)- A/" (X r - x) l t < T „] 

= ±E [A/' 1 (X Ta - x) l t>TR ] > on (0, oo) x E. 

that is u is a supersolution of the heat equation. Further, equation (17) applied with t = oo and using that 
X TR ~ i/ also gives 

u n (t, x)-u v {x) = ^j E [A/™ (Xr - x) l r < TR ] dr. 
Now for (t, x) £ R, also (s, a;) € R for s > f and one expects that A/" (XJ R — x) w a.s. since for n big enough 

suH»{A/ n (--x)}w{x}, 
hence the Ito correction vanishes and above reads in the limit as n — > oo as 

u = u u on R. 

Now if (t, x) £ R c and t > tr then the properties of the Root barrier imply that X TR ^ x, hence A/ n (X TR — x) l t> 
for n big enough. We therefore arrive (as n — >• oo) at 

d t u = Au for (t,x) G i? c , 
u(0, x) = u M (i). 

To conclude this (very formal) argument note that RlJ R c = [0, oo] x M and therefore putting all the above 
shows us that 

min (u — u v , dtu — Au) = on (0, oo) x R, 

u(o,-) = 

The proof of Theorem 5 makes this type of reasoning precise via the method of semi-relaxed limits. 

Theorem 5. Let (a,fi,v) be either as in Theorem 3 or as in Theorem 4- Let R G 1Z (a, /it, V) and denote with 
u the potential function of the stopped process X TR = (X t ^ TR ) t , i-e. 

u (t, x) := u x -r (t, x) := -E [\XJ R - x\] . 

Then {u} = Voo (c, Ufi, u v ), in other words, the potential function of X TR is the unique viscosity solution of 
linear growth of the obstacle problem (10) with initial data u^, barrier u v and heat coefficient a 2 , 

,^g, ( min (u — u u , dtu — ^-Atij = 0, on (0, oo) xM 

\ u(0,.) = u M (.). 

Moreover, 

(1) for every xglfHM (t,x) is non- increasing and 

u u (x) < u (t, x) < u M (x) V (t, x) , 

(2) u\ r =u v \r andV(t,x) lim^oo u (t, x) =u v (x), 

(3) u is Lipschitz in space (uniformly in time), 

\u(t,x) -u(t,y)\ 
sup sup j j < oo. 

te[o,t») x^ y \x — y\ 

Proof. We have to show that u is a viscosity solution of (18) and to show uniqueness of the solutions via the 
comparison theorem given in the appendix we also need to establish that u is of linear growth. First note, that 
to verify that u is the viscosity of (18) it is enough to show that (in viscosity sense) 

(19) 



u — 




> 





on 


[0,oo) x 


^. 












)» 


> 





on 


[0, oo) x 


u — 









on 


R. 














). 







on 


R c , 
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IX, 



sc|] = u v (x) 



The first inequality follows immediately via conditional Jensen 

(20) u (t, x) = -E [\X T t R - x\] > E [E [- \X TR - x\ \Ft hrR }] = 

For the other inequalities we approximate (18) and conclude by the Barles-Perthame method of semi-relaxed 
limits. 

Step 1. u n (t,x) := E [-0™ (X T t R - x)\. Define the sequence (ip n ) C C 2 (1 

r x ry 

ipn (x) = / / ncj) (nz) dzdyMx £ M. 
Jo Jo 

where d> is the usual Gaussian scaled to the unit disc 



as 



4>{x) 



exp 




l-x 2 



for \x\ < 1 
otherwise. 



Especially note that tp n (.) — > |.| uniformly, Aij) n (.) is continuous, < Aijj n < O (n) and supp (Aip n ) C \—h> 1 
(we could replace i\) n by any other sequence with this properties). Further introduce 

u n (t,x) = -E [tp n (Xl R - x)] , 

K (x) 



Since ip n (.) — > |.| uniformly we have P-a.s. 

\r (^ T "--)-K rfi --IL ;[ o,oo)x: 
hence we get uniform convergence of u n ,u™ and u™, i.e. 
(21) l« n -«U[o,=o)x 

\K- u Aoo, 



i> n {y-x)ii (dy) . 



sup \r (XJ« -x)- \Xl 

(t,a:)e[0,oo)xR 



X —> r . 



Denote the complement of R with 

R c = ([0,oo] x [-00,00]) \i? 
and note that since R is closed, R c is open. Further, by dominated convergence we have Vx £ R, Vn 6 N that 
(22) lim u (i, x) = Un (x) and lim u n (t, x) = u n (x) . 

t— foo t— »oo 

Step 2. u - = on J?. Fix ifl and apply the Ito-formula to — n (. — x) and the local martingale X TR . 
After taking expectations and using Fubini's theorem we arrive at 



(23) 



u n (t, x) 



o 2 (r,X r ) 



A^ n (X r - x) l r<7 



dr. 



Taking lim^oo on both sides in the above implies by (22) that 

V 2 (r,X r ) 



(24) 



ul (x) = ul (x) 



E 



Aip n (X r - x) l r<7 



dr. 



Since ^-A^™ > 0, (23) shows also that u n is non-increasing in time and by subtracting (24) from (23) we arrive 
at 

V 2 (r,X r ) 



u n {t,x)-ul (x) 



2. 



-Alj) n (X r - x) 1 



r<r B 



dr > 0. 



The uniform convergence given in (21) implies that t i— >■ u(t, x) is non-increasing, i.e. \/(t,x) £ [0, oo) x R we 
have that 

u(t + r, x) < u(t,x) Vr > 0. 

Fix a point (t,x) £ R° (R° denoting the interior of R). Because R is a Root barrier this implies that Vr > t 
also (r, x) £ R° and since R° is open we conclude that 3uq s.t. 

1 



(r, y) £ [0, oo] x [— oo, oo] : r > t, \x — y\ < 



C R°. 



Since supp (Aip n ) C [— n 1 ,n - 1 ] this shows Atp n (X r — x) l r < TR = for all r > t and all n > no, hence 

V 2 (r,X r 



lim (u n -u™)(t,x) = 



lim A0 n (X r - x) l r 



dr = V(t,x)c R°. 
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By continuity of u and u v we can conclude that 



u\r= u,/ In- 



Step 3. (dt — ^-AJ u > on [0, oo) x E. From (23) it follows that u n (t, x) has a right- and left- derivative 



V (t, x) £ [0, oo) x R; to see this take 



d t +u n (t,x) = lime 1 (u (t + e, x) — u (t, x)) 



-E 



lime 

<J 2 (t,x t ) 



^ (r,X r ) 



Aijj n (X r - x) l r <r R dr 



AV>" (X t - x) l t 



similarly it follows that the left derivative dt~u n is given as 



(25) 



d t -u n {t,x) 



-E 



a 2 (t,X t 



■A^ n (X t - x) l t 



negative we conclude 



Note that for every t <E [0, oo) d t -u n (t, .) , d t+ u n (t, .) G C°° (E, E); further, since ^ ^ Xt) Aip n (X t - x) is non- 

dt-u 11 (t,x) < d t+ u n (t,x) < 0. 
(Note that for n < oo we cannot expect the above inequality to be an equality; e.g. consider v = <$o, /i = Af (0, 1), 
(7 = 1 which is solved by R = {(t,x) : t > 1}, hence tr = 1). From the definition of u n it follows that we can 
exchange differentiation in space and expectation to arrive at 



(26) 



° (*' ^ Au n (t, x) = - ° 2 — ^ E [A-0" (X t r « - x)] < on [0, oo) x E, 



2 v ' ' 2 

which is a continuous function in (t, x). Lemma 5 shows that V(a,p, m) <E V'q U 11 (t,x 

a > d t -u n (t, x) and m < Au n (t, x) . 

Hence, by (26) and (25) 



a — m > dt-u n (t, x) 
1 



cr 2 (i,x) 



Aw n (t, x) 



E <H (i, x) A^ n [XI* -x)-a z (i, X t ) AiP n (X t - x) l t 



<tr\ 



Splitting the term inside the expectation gives 

dt^u t (t 1 x)~m > — E [cr 2 (t, x) Aip" (X£ R — x) — a 2 (t, x) Aip n (X t — x) l t < T 
2 v .. 



+i E [a 2 (i, x) A^ n (X t - a:) l t < TR - a 2 (t, X t ) A^" (X t - x) l t < TB ] 
- >> .. ' 



=:II n (t,x) 



-I n (t,x) + -II n (t,x). 



We conclude that u n is a supersolution of (j) t — irA^j u — ^ (/„ + iT„) = on (0, oo) x E. Further, 

I n (t,x) = ^[a 2 (t,x)Aip n (X^ -x)-a 2 {t,x)Aip n (X t -x)lt< TR ] 
= cr 2 (t, x) E [Atp n {X TR - x) l t>TR ] > 
hence u n is also a viscosity supersolution of 

d t u - slAu - \ll n = on (0, oo) x E 

tt(o,.) = «»(.). 

Using the Lipschitz property of cr, supp (Aip n ) = [— and that |A^> n | 0O < cn we estimate 



\II n (t,x)\ < E[\a 2 (t,x)-a 2 (t,X t )\ Axp n (X t - x) l t 



<tr\ 



< VLip — Z&LG 

n 

< ^Lip^LG 



l + x + - E [Aip n (X t - x) l t<TR ] 
n 

-E [cnl\ Xt -x\<n-i] 
¥(\X t -x\ < n- 1 ) 
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(<JLip and <jlg as defined in (4) and for the second inequality we use the trivial estimate | cr 2 (t, x) — a 2 (t, X t )\ < 
|<7 (t, x) — a (t, X t )\ \a (t, x) + a (t, X t )\ combined with \x — X t \ < — , Lipschitzness and linear growth of a). Now 
for every compact K C [0, oo) x R we have 

lim sup P \\X t ~x\< n,- 1 ] = 
n ^°°(t,x)eK 

since by Assumption 1 the process X has a continuous density with respect to the Lebesgue measure on (R, B) 
and 

r x +i: 

PO-Xt-sl^n- 1 ] = / \y-x\f(t,y)dy 



< 



— I f (t, y) dy — > as n — > oo 



uniformly in (t, x), therefore //„ —> locally uniformly on [0, oo) x R. By step 1, u n —> u and u n (0, .) — > (.) 
locally uniformly as n — > oo. The stability of viscosity solutions, Proposition 4, implies that u is a viscosity 
supersolution of 



(27) 



r (ft-^A)ti 

1 "(0,.) 



u = on [0, oo) x 

•) = « M (0 
which proves the desired inequality. 

Step 4. ^<9t — ifAj u = on i? c . From Lemma 5 it follows that if dt-u n (t,x) < dt+u n (t,x) then 
p2,+ u n = ( m which case we are done) and if dtu n (t,x) = dt~u n (t,x) = dt+u n (t,x) then V(a,p, m) G 
p2,+ u n (t^ x ), a = d t u n (t,x) and m > Au n (t,x). Hence, in the latter case we have V(a,p, m) E V 2 ' + u n (t,x) 
that 

a-m < d t u n {t,x)- a ^ ,X ^ Au n {t,x). 
Proceeding as in step 3 we see that u n is a subsolution of 

d t u - ^Au - | (7„ + //„) = on (0, oo) x R 



u(0,.) = «»(.)• 

We already know that //„ — > locally uniformly and now show that I n — > locally uniformly on R c : since R 
is a Root barrier we have 

{T R + r,X TR ) E R Vr > 0, 
hence if (£, x) E i? c and t > tr then X TH 7^ x. Therefore 

lim sup Aip n (X TR — x) 1{>t h = for every compact K C R c 

which is enough to conclude that I n converges locally uniformly on R c to 0, i.e. for every compact K C R c 



lim sup I n (t,x) < \<t 2 \ limE 



sup Aip n (X TR - x) l t >r R 

(t,x)EK 



= 0. 

The stability results, Proposition 4, implies that u is a subsolution of 

d t - y A^j u = on R c , 
«(0,-) = «„(•)• 

Together with the result from step 3 this shows that u is a solution in viscosity sense. 

Step 5. Putting the above together shows that u is a viscosity solution of the obstacle problem (18). 
To see that u is of linear growth, recall that by definition u(0,x) = u^(x) and in step 2 we showed that 
Vx E R, t 1— ► u (£, x) is non- increasing. Further, we know u (t, x) > (x) V (i, x) from the estimate (20), hence 
\u (t, x) < I u M (x)| + \u v (x)|. This together, with the properties of potential functions, Proposition 1 and Lemma 
2, implies sup( t x ) e [o oo) x k i+|if < 00 • This shows that the assumptions of the comparison result, Theorem 9, 
are met and we conclude that {u} = Voo (c, u M , u v ). □ 

For a function v E C 1,2 ((0, 00) x R, R) one has P 2 ^(f,x)nP 2 ^(M) = { (ff !(*,*), ^| (*,»). &l(^)) }■ 
The lemma below is a generalization of this for functions which are only left- and right-differentiable which 
appear in the proof above. 
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Lemma 5. Let v G C ((0,oo) x R, R) and assume that V(i,x) G (0, oo) x R, v has a left- and right- derivative, 
i.e. the following limits exist 

v(t + e,x) -v(t,x) v(t + e,x) - v(t,x) 

a t +v (t, x) = hm and o t -V (t, x) = lim , 

e\o e eyo e 

If dt-V (t, x) < dt+v (t, x) then 

a G [d t -v (t, x) , d t +v (t, x)] V (a,p, m) G V 2, ~v (t, x) . 

Ifdt-v(t,x) <dt+v{t,x) then V 2 ' + v (t, x) = and if dt~v (t,x) = dt+v(t,x) then V (a,p,m) G V 2 ' + v (t, x), 
a = d t v(t,x) (= dt~v(t,x) = dt+v(t,x)). In all the above cases, if v is additionally twice continuously 
differentiable in space then 

m < Av(t,x) V (a,p, m) G V 2 ~v (t, x) , 
m > Av(t,x) V(a,p,m) eV 2 ' + v(t,x) . 

Proof. Every element (a,p,m) G V 2 '~v (t, x) fulfills 

v (t + e, x) — v (t, x) > ae + o (e) e — > 0. 

Applied with a sequence e n /* 0, it follows after dividing by e™ and letting n — > oo that dt~v(t,x) < a. If 
9t_w (i, x) < 9t+w (t, x) then for (o, p, m) G P 2 ^?; (i, x) we have 

v (t + e, x) — v (t, x) < ae + o (e) 6 — >• 

which leads after taking e'/O resp. e" \ to 

dt+v (t, x) < a < dt-V (t, x) 

and hence contradicts the assumption dt-v (t, x) < (t, x). The other statements follow similarly. □ 

Remark 10. The proof of Theorem 5 holds under weaker assumptions, e.g. only assumption 1. 

4.2. Prom the obstacle problem to the Root barrier. Recall Theorem 3 and Theorem 4 which ensure 
that if the measures are in convex order, a fulfills linear growth and Lipschitz conditions and 1 1— ► [X] t increases 
fast enough, then there exists a Root barrier, i.e. TZ (a, (J,, v) ^ 0; and further, that under the conditions of 
(/x, ^-regularity there is a unique barrier, i.e. \R- reg (a, /i, v)\ = 1. The theorem below shows that the unique 
(/i, ^)-regular Root barrier can be calculated by solving the obstacle PDE (18). 

Theorem 6. Let (a, /i, u) be either as in Theorem 3 or as in Theorem 4- Then there exist a unique viscosity 
solution of the obstacle problem (18), i.e. {u} = Voo (<r, u^, u u ), and 

(28) {(t, x) G [0, oo] x [— oo, oo] : u (t, x) = u v (x)} = R 

where {R} = lZ reg (a, fi, v), i.e. R is the unique ([A, v) -regular Roob barrier which solves the Skorohod embedding 
problem given by (o~,[i,v). 

Proof. Above assumptions guarantee that [a, fx, v) fulfills the conditions of Proposition 5, i.e. there exists a vis- 
cosity solution of the obstacle problem with initial data u^, barrier u v and heat coefficient a 2 , i.e. Vt (c, U^, u v ) = 
{ u I[o,t]xk} f° r every T G (0, oo). On the other hand Theorem 3 gives the existence of a unique (/x, i/)-regular 
Root barrier, i.e. \lZ reg {a, fx, v)\ = 1 and from Theorem 5 we know that 

(t,x) i y -E[\x - X tATR \] 

is the unique element of Voo (c, u^, u u ), i.e. is a viscosity with linear growth uniformly in time. This is sufficient 
to conclude that u (t, x) = — E [|x — A tATR |] for all (t, x) G [0, oo] x R, i.e. we have shown that 

{U} = Voo (cTjMpjWy) . 

It remains to show (28). From step 2 in the proof of Theorem 5 it follows that 

R C {(t, x) G [0, oo] x [— oo, oo] : u (t, x) = u v (x)} =: Q. 

To see the other inclusion Q C R note that by definition of u we can apply the Tanaka formula to see 

u (t, x) = Ufj, (x) - E 

where if denotes the local time of the stopped process X T = (A TAt ) 4 . Hence 

u (t + r,x)-u (t, x) = E [Lj — L* +r ] . 

Now if (t, x) G Q, then (t + r,x) G Q and so u (t + r, x) — u (t, x) = u v — u v (this follows from properties of the 
PDE solution as stated in Theorem 5). Hence 

LI = L* +r V(t,x)GQ. 

Therefore it is sufficient to show that if (t,x) G R c then 3r > s.t. 

P [L* +r > Lf ] > 
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However, as a consequence of the Ray-Knight theorem [24, Chapter 22, Corollary 22.18] we have V (t, x), r > 
that 

{L* +r - L* > 0} = J inf XJ<x< sup Xj\ 

[se[t,t+r] se[t,t+r] J 

so it only remains to to show that V (t, x) £ R c there exists a r > s.t. 

inf X J < x < sup AJ > 0. 

If the assumptions of Theorem 3 hold, i.e. X a geometric BM this follows immediately and if we use the 
assumptions of Theorem 3 we can use Dambins/Dubins-Schwarz to see that X = Xq + W\x] for some BM W 
on (fi, J 7 , (J" t ) , P), and since P [tr > 0] > we have 



inf Xq + W[ X ] At r < X < SUp Xq + W[ X ] Ar fl 
»e[i,t+r] sf z[ t .t+r] 



> 0. 



□ 



Remark 11. Let Q € 1Z (cr, fi 7 v). Applying Theorem 5 with Q leads to a function u (t, x) := E [| X T t Q — x|] , u € 
V (cr, w^, u„) and applying Theorem 6 to it yields the barrier QuN^' v . In other words, the PDE representation 
always chooses the unique (/i, i^)-regular barrier. 

Remark 12. Consider (cr,/i,^) = (l t <i,(5o,A/'(0, 1)) and note that a(t,x) = l t <i does not fall into our setting. 
Obviously any fixed (deterministic) r > 1 is an element of S (cr, /i, v) but only r = 1 is the minimal solution to 
the embedding problem (in the sense of Definition 6). Note that the above construction via the obstacle PDE 
chooses the minimal solution r = 1. 

4.3. The reversed Root barrier (Rost's barrier). Root's barrier R resp. the associated hitting time lets 
the process X diffuse as much as possible before it hits the barrier. Motivated by this, one can ask for the 
other extreme: how to build a barrier R such that its hitting time solves the Skorokhod embedding problem 
and stops the process X as soon as possible. Indeed such a construction was given by Rost [32] who was led 
to this topic by the Chacon-Ornstein ergodic theorem and an extension of the filling scheme from discrete to 
continuous time. A closed subset R of [0, oo] x [— oo, oo] is a Rost barrier R if 

(1) (t, x) e ~R implies (s, x) E ~R Vs < t, 

(2) (0,x) 6 ~R Vx £ [-00,00], 

(3) (i,±oo) 6 R Vt E [0,+oo]. 

Following the very same approach as in Theorem 5 via semirelaxed limits allows one to show that if R is a Rost 
barrier which solves the Skorokhod embedding problem for (cr, fi, v) then the continuous function 

u (i, x) := u (t, x) := ux t ~r. (t, x ) — u v {x) 
is the unique viscosity solution of the PDE 

f d t u = min (0, 4 Aw ) » 011 (°> 00) x E 
I u (0, •) = (•) - u u (•) . 

However, while the PDE gets simpler — in the sense that standard parabolic comparison results apply and 
existence follows from Perrons' method — the question which conditions on (cr, /1, v) actually guarantee the 
existence of a Rost barrier which solves the Skorokhod embedding problem and which additional assumptions 
on the Rost barrier guarantee its uniqueness is much less clear to us; i.e. we are not able give an analogue of 
Theorem 6 which applies to general (cr, fx, v). In fact, the existence of the reversed/Rost barrier has been studied 
only for certain cases by Rost [34] (as well as [31, 32, 33] and we point to forthcoming work 5 of Cox et al. for 
more sufficient conditions). 

5. Application i: calculating the root barrier via the barles-souganidis method 

Theorem 6 allows to calculate the Root solution of the Skorokhod embedding problem, {R} = 7Z reg (cr, fi, v), 
by calculating the viscosity solution of the obstacle problem (10) with initial data u^, barrier u v and heat 
coefficient cr 2 . One of the benefits of the viscosity approach is that one can resort to the Barles-Souganids 
method [3, 4] when discussing numerical schemes. (The Feynman-Kac representation of the solution to (2) 
as the solution to the RFBSDE (13) serves a secondary purpose, namely, that one can make use of existing 
numerical methods for the RFBSDE to obtain an approximation of the Root barrier). In fact, due its relevance 
for optimal stopping problems there exists a vast literature on schemes for the obstacle PDE (18). While it falls 
outside the scope of this article to study numerics of the obstacle PDE (10) in full generality we still give two 



(29) 



5 Personal communication. 
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brief applications: firstly we show that the classic Barles-Souganidis convergence result can be easily adapted 
to give convergence when /x, v have compact support, secondly, we discuss the case of general [i, v but under the 
extra assumption a = 1 (i.e. embedding into Brownian motion) in which the literature on viscosity solutions 
even provides good rates of convergence. 

5.1. [i and v of bounded support. The assumption of compact support implies (via Proposition 2) that 
(N fJ " l/ ) c C [a, b] for some a, b G R, hence (10) reduces to a PDE on a compact set with uniformly continuous 
<r, u^, u u data. Under the assumptions of Theorems 3 or 4 we have {u} = Vqo (a, u^, u v ) and the results in [3, 4, 1] 
(see also [17]) can be easily adapted to prove convergence of numerical schemes. We give a quick construction 
using an explicit finite differences scheme to approximate the solution of (2). On Ot '■= [0, T] x [a, b] and setting 
h := (At, Ax) = (j^;, j^r^J for Nt, N x G N large enough we define the time-space mesh of points 

G h := {t n : t n = nAt, n = 0, 1, • • ■ , N T } x {xj : x 3 ■ = a + jAx, j = 0, 1, • • • , N x } . 

Let B(O t ,R) be the set of bounded functions from O t to K and BUC(O t ,R) C B(O t ,R) the subset of 
bounded uniformly continuous functions. Take -0 G BUC (0r,R), we define its projection on Qh by ip h : O t — > 
R with O t := [0,T + At) x [a - Ax/2, b + Ax/2] as ip h (t,x) := ip(t n ,Xj) when (t,x) G [t n ,t n+1 ) x [xj - 
Ax/2, Xj + Ax/2) for some n G {0, 1, ■ • ■ , Nt} and j G {0, 1, • ■ • , N x }; of course i/j h G B [O t , K). We denote 
the approximation to our solution u G BUC{O t ,M) of (2) by u h G S(O t ,K). Below we use a standard finite 
difference scheme, i.e. we define the operator S h : B {Ot, R) x [0, T] x [a, 6] i— y M. as 

f U/( (a?) ^ ,(t,x) G [0,At) x (a,b) 

S h [u h ] (t, x) := I u h (t, x) + A 2ff ) f ) (u h (t, a; + Ax) - 2u h (t, x) + u h (t, x - Ax)) , (t, x) G [At, T] x (a, b) 
{ Up{x) ,{t,x) G [0,T] x {a, 6} 

where the well known CFL condition (computed for this case) CFL := l^loo-Ta. 6]x[o r] < (Ax) 2 is assumed 
to hold. The values of u h are computed by solving for u h (t,x) via G(.) = where G : (0,oo) 2 x Ot x M x 
B (O t , K) h> R is defined as 

G(/i, (f + Ai, .t) , u^' (t + At, x) , m' 1 ) := min{u' 1 (t + At,x) - u v (x) , u h (t + At, x) - S h [u h ] (t,x)} 

(for (t + At, x) G Ot of course). Following [ 3, 4] we only have to guarantee that the operator S h [.] (.) defined 
above and the PDE satisfies along some sequence h := (At, Ax) converging to zero the following properties 

• Monotonicity. G (h,(t,x) ,r, f h ) < G (h,(t,x) ,r,g h ^ whenever / < g with f,g G B (and for finite 
values of h, t, x, r); 

• Stability. For every h > 0, the scheme has a solution u h on Qh that is uniformly bounded independently 
of h (under the CFL condition, see above); 

• Consistency. For any ip G C£° (Ct;R) and (t,x) G Ot, we have (under the CFL condition, see above): 

hm ((^ (i„ + Af,X;,-) + - A — 



^ (t, x) - u„ (x) , f StV - yA^J (t, x) 

• Strong uniqueness, if the locally bounded USC (resp. LSC) function u (resp. v) is a viscosity subsolution 
(resp. supersolution) of (2) then u < v in Ot; 

Proposition 6. Let T G (0,oo). Assume \i,v have compact support and that o~,\i,v fulfill the assumptions of 
Theorem 6. Then u h G B([0,T] x R, R) and 

K - U L;[0,T]xR^° ash -KM) 

where u is the unique viscosity solution of (2) on [0,T], {it} G Vt (cr,/x, f). 

Proof. This results is more or less a direct consequence of [3, Theorem 2.1] and [1, Theorem 2.1]: the required 
strong uniqueness of the obstacle PDE (2) follows from the comparison theorem in the appendix while the exis- 
tence from Theorem 5. The monotonicity, stability and consistency follow (stability and consistency conditions 
of the scheme hold under the CFL condition above) by a direct calculation (which we leave to the reader; see 
e.g. [4, 3, 1, 2]). The rest of the proof is a direct adaptation of the proofs [3, Theorem 2.1] and [1, Theorem 2.1] 
combined with the remarks on the first example in [I, Section 5]. One first shows that the operator S h [.] (.) 
approximates indeed the PDE component of (2) then one adds in the barrier to recover the full equation (2). 
The core ingredients are the semi-relaxed limit technique of Proposition 4 combined with the strong uniqueness 
result (see Theorem 9) as sketched out in [1, pl30]. □ 
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K(t, X) 




l.Of 




Figure 1. a = 1, /i = 6 and v = ±6_i + §<5 + |<5i (hence iV^" = M \ (-1, 1) U {±oo}). 
Applying above explicit finite difference scheme with CFL = 0.2 and 50K time steps on the 
interval [0,2] gives an approximation of (2) on the spatial domain [—1,1]. The figure on the 
right shows that for to = 0.39348 we have u (to, 0) = u v (0) = | which determines the spike at 
x = for the Root barrier depicted on the right. 





FIGURE 2. a = = S and v = + \5_l + |§<53. 
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FIGURE 3. a (x) = x, fi = S\ and v = U |1) . Note that while target and initial distribution 
are symmetric, the barrier becomes asymmetric due to a (x) = x. 




FIGURE 4. a (x) = 1,/x = f 5_i + \5 .5 and v = i<5_i + |W ([0, 1]) . 



Remark 13. An issue is the possibility of smoothing and u v prior to the numerical approximation. Usually 
u^Un have "kinks", for instance a common setup is with [i = So leading to u^x) ~ — \x\ with a kink at x = 
or v has atoms and these atoms are (as expected) a source of instability for finite difference schemes. Within 
the financial framework, the measure v is given by market data, the observed prices of call options, via the 
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Breeden-Litzenberger formula. In practice the data available is only for finite amount of strikes and hence the 
inferred measure v will have atoms. It seems intuitive that smoothing u v (or v) leads to a more stable method 
(see for example [19]). What is not clear is how to smooth u v (or v) without adding arbitrage. 

5.2. Calculating general embeddings for Brownian motion (/i and v are of unbounded support). 

We now discuss the case of general [i, v. For simplicity we restrict ourselve to embeddings into Brownian motion 
(i.e. a = 1) in which case more recent results of Jakobsen [22] directly apply and lead to a rate of order 1 (if 
one is only interested in convergence the argument given above can be modified). Denote again h = (At, Ax) 
and consider schemes of the type 

Uh (t + At, x) = max (S&tUh (t, x) , u v (x)) 

where SAt is the (formal) solution operator associated to the heat equation dw — |Aw = 0. In the case that 

we use a finite difference method to approximate S the scheme can be written as 

(30) 

( Uh(t + At,x) - u h {t,x) u h (t, x + Ax) - 2u h (t, x) + u (t, x - Ax) . . \ 

mm — ——2 , Uh (f + At, x) - u v (x) = 0. 

V At ( Ax ) J 

A direct calculation shows that this is equivalent to 



Uh (t + At, x) = max 




This representation can then be used to prove the following 

Proposition 7. Let ji < cx v, (i,v € M 2 m - There exists a unique Uh solving (30). Further, if At < (Ax) 2 , u 
denotes the unique viscosity solution of (2) ({u} = Vt (1,?V, ■«„),) and uo-h is an approximation of uq which is 
bounded independently of h then 

\ u - u h\oc~ SU P \ u ~ u o-,h\ + [Ax) 1/2 

[0,A*)xR 

Proof. This is a direct consequence of [22, Section 3]. The main idea is to show that one can replace the 
Barles-Souganidis assumptions by more special conditions (CI — C5 in [22, Section 2]) which can be checked 
for the finite-difference scheme. □ 

Remark 14. In any real life implementation one is forced to deal with a domain truncation procedure in order to 
approximate the solution, the convergence of the truncated PDEs has to be shown as well as existence/uniqueness 
of the solution to the problem in the space of discontinuous (due to truncation) viscosity solutions. Moreover, 
the truncation creates the appearance of an artificial boundary at the truncation points and in this boundary one 
can impose either Cauchy- or Neumann-type conditions for the approximation leading to non-trivial questions. 
Note that under the compact support assumption for fj., v treated in the first part even though one interprets 
the PDE in a bounded domain, u^ and u v are known outside this domain hinting that one could use higher 
order finite differences stencils for the spatial derivatives without having to introduce an artificial boundary as 
in the usual case. 



6. Application ii: model-independent bounds on variance options 

We briefly recall an application of Root's barrier to obtain model-independent bounds on variance options 
(which was Dupire's [ original motivation for looking at Root's solution): much of financial mathematics is 
concerned with the pricing of claims contingent on the evolution of a risky asset. One approach is to 

• firstly, postulate a model for the risky asset, viz., denote with (Q, J-, (J- t ) ,P) a probability space satis- 
fying the usual conditions carrying a real- valued process P = (P t ) which describes the evolution of the 
price of a liquid asset and with constant interest rates r the forward price process S = (St) of this risky 
asset is given as 

St = e~ rt P t - 

• secondly, calculate an arbitrage-free price in this model as the expected value of the claim under a 
risk-neutral measure, viz., if Ft = F (St : < t < T) denotes the payoff of the claim, one considers a 
(so-called pricing) measure Q equivalent to P under which S is a local martingale and an arbitrage-free 
price is then given as E 1 ® [Ft] ■ 

Indeed, if Ft equals (St — K) + or (K — St) + these claims form a liquid market and from observed prices one 
can in principle calibrate the chosen model. However, this approach falls short when Ft is an exotic, non-liquid 
option and in practice a trader is faced with the following question: how to price a given 



F T = F(S t : < t <T) 
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(which might be highly path-dependent hence prices might be highly model-dependent) under the only as- 
sumptions that S is a local martingale and that prices c (T, K) for Ft = (St — K) + are known for every T, K. 
Surprisingly there is a straightforward link with the Skorokhod embedding problem — see the excellent surveys 
of Obloj [28] and Hobson [21] — and we recall the argument in a nutshell: 

• the market prices c (K, T) give (via the Breeden-Litzenberger formula) the distribution of the increment 
of the forward price, i.e 

fc(J\x) 

Ltiax) = r-= — ax, 

ax z 

• any calibrated model (Q, F, (Ft) , P) with forward price process S must have St ~ So ~ /i, 

• assuming Dambis/Dubins— Schwarz applies we can write S, — So = W T r\ (W denoting a Brownian 
motion, r (.) a time-change), hence 

W T(T ) ~ A 1 

and r (T) solves the Skorokhod embedding problem 6 . 
From the above we conclude that all arbitrage-free prices of Ft are contained in 

inf F (W T ( t ) :0<t< T) , sup F (W T ( t ) : < t < T) 
r r 

with sup and inf taken over all stopping times r which solve the Skorokhod embedding problem for [i, further, 
the bounds are attained by "extremal" solutions of the Skorokhod problem which in the case of variance options 
is given by the Root and Rost solutions. 

We now give a rigorous argument for the latter case of variance options: henceforth assume without loss of 
generality that So = Po = 1 and that 

F T = f ([In P] T )=f ([In S] T ) 

where / is a convex, increasing function (e.g. / (.) = (. — k) + for fixed k > 0). 

The interest in the Root model is that it gives the universal lower bound among all possible arbitrage-free 
forward price processes S with a prescribed marginals. 

Theorem 7. Let S be a positive, real-valued continuous local martingale carried on some probability space 
(fl,F, (Ft) , Q) satisfying the usual conditions. Further assume, [InS]^ = oo a.s. and 

St ~ v. 

Let R G TZ (id, S\, v) . Then for every f : [0,oo) — > [0,oo), / (0) = which is convex and has a bounded right 
derivative, 

(31) E[f(r R )]<E®[f([]nS} T )]. 

Moreover, this bound is sharp in the sense that there exists a continuous local martingale S R (which is indepen- 
dent of f) such that (31) becomes an equality. 

Proof. It follows from Lemma 7 given below that there exists a Brownian motion B on (0, T, (Ft) ,Q) such 
that 

[in sy 



St = cxp ^-B[inS] t - 

Since by assumption St ~ v, and (£t/\[inS] T ) t i s u -i-; it follows that [kiiS] T is a solution of the Skorokhod 
embedding problem S (id,8\,v). Now using the minimality of the Root solution, Theorem 4, implies (31). To 
see that this bound can be attained consider a filtered probability space (f2, F, (Ft) , P) satisfying the usual 
conditions which carries a Brownian motion B. Using the Root barrier R we define a process on this space 

St ■= cxp \B T ,t - - ( tr A 



t r a T '_ t 2 \ ±h T* t 

Then S R is a continuous, local martingale with respect to the filtration JF R := J 7 t . Sj$ ~ v and 

[lnS R ] T = T R . 

That it, S R is an admissible pricing model which matches the marginals and [lnS' ii ] T attains the lower bound. 



□ 



Lemma 6. A continuous, positive process S on (Q,F, (Ft) ,Q) is a continuous local martingale (fl,F, (Ft) 
iff 

S. = So cxp (m. - i [M]) a.s. 
for some continuous, local martingale M . In that case M. = j S~ 1 dS a.s. 



"To make this step rigorous needs a bit more care since there might not exists such a time change etc. 
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Proof. Follows from a direct application of the Ito formula to In St ■ □ 

Lemma 7. Let S be a positive, continuous local martingale on (Q, T, (Tt) , Q). Further assume pnS] = oo 
a.s. Then there exists a Brownian motion B on (£l,T, (Tt) , Q) such that 

[In S]' 



S. = S cxp I B [ i nS] 

Proof. Let M be the continuous local martingale M of Lemma 6. Since In S. = In So + M. — | [M] we have 
[-^Hoo = [ m ^]oo = 00 a - s - By ^ ne Dambis/Dubins-Schwarz Theorem (cf. [29, Chapter 5, Theorem 1.7]) there 
exists a Brownian motion B on (SI, T, (Tt) , Q) s.t. M t = -B[m] Vt, Q- a.s. □ 

Remark 15. The assumption [InS 1 ] = oo is needed for the Dambis/Dubins-Schwarz Theorem to guarantee 
the existence of a Brownian motion B on (fl,T, (Tt) , Q). It can be dispensed with for the price of using an 
enlargement of (0, T, (T t ) ,Q) cf. [29, Chapter 5, Theorem 1.7]. 

Remark 16. Above proof does not easily generalize to match more marginals for general payoffs /: for example 
consider the two-period case (to,t\,t2) = (0,1,2) with given marginals (S\, /xi, /X2). For every given local 
martingale 5, [ln5] 1 <E S (id, 5i, /Xi) and [ln5] 2 — [ m <51i G <5 (id, 1^2). By the superaddivity of / one gets 
a lower bound E[/ ([ln5] 2 )] > E[/([lnS] 2 - [lnSy] + E [/ ([In S}^], however the right hand side can be 0. 
For example, let fi± and fj,2 be the distribution of a geometric Brownian motion at time 1 and 2, in which 
case [ln5] 1 = [hi5] 2 — [lnS'] 1 = 1 and if one takes /(.) = (. — 1) + the above reduces to the trivial estimate 
E[/([lnS] 2 )] = 2>0 = /(l) + /(l). 

7. Appendix: existence and uniqueness of the root solution for diffusion martingales 

As pointed out in Subsection 2.3, the results in this section are generalizations of results which cover the case 
(a,p,v) = (l,<Ji,/x), cf. [30, 26, 21, 28]. However, we could not find precise statements in the literature which 
apply to the case of general (a, (J,, v) with /x, v of possibly unbounded support' and since for example some issues 
arise when a is of linear growth (resulting in the possibility of [X] < 00 a.s. which causes some difficulties) 
we decided to give a complete argument for the convenience of the reader. 

We first recall a couple of results on Root barriers and the relation with their first hitting times. The first 
result concerns the metric used in the barrier space. 

Lemma 8. Equip [0,1] x [1,1] with Euclidean metric and denote with r the metric on [0,+oo] x [— 00, +00] 
induced by the homeomorphism 

[0, +00] x [-00, +00] -> [0,1] x [1,1] 



(t,x) 1 y 



t 



1 + t 1 



Then, 

(1) the map 



p (R, Q) := sup r (x, R) V sup r (x, Q) 



is a metric on the space C of closed subsets of [0, +00] x [—00, +00], 
(2) (C,p) is a separable, compact metric space and 1Z is a closed (hence compact) subset. 

Proof. Root [30] (see also Loynes [26]). □ 

With this metric, the map that associates a barrier to its hitting stopping time is uniformly continuous from 
(TZ, r) into the set (T, c.p.) of first hitting to barriers equipped with the convergence in probability (c.p.). The 
next lemma states this more precisely. 

Lemma 9. Let Assumption 1 hold. Let B £ TZ with corresponding first hitting time Tb- Iff [t~b < 00] = 1, then 
Ve>03<5>0 s.th. if R £ TZ with tr its corresponding hitting time and r (B, R) < S then P [tr > tb + e] < e. 

Moreover, take a sequence of barriers (R n ) n an d the associated hitting times (r„) „, if the sequence of bar- 
riers converges (w.r.t. the metric r), R n — >„ R and (defining tr is the hitting time associated to R and) if 
sup n E [r n ] = K < oo, then it holds that E [tr] < K and P [|tr„ — tr\ > e] — > . 

Proof. The proof is a direct reformulation of the proofs of Lemma 2.4 of [30] and Lemma 1 of [26] which were 
stated with the Brownian motion being the underlying time-space process. Just remark that under Assumption 
1 the solution of (6) has continuous paths, the strong Markov property, no interval of constancy and covers the 
required state space. □ 

We can now start proving the main theorems. 



'An excellent article is Rost's [34] but he just treats Markov processes taking values in a compact space. 
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7.1. Existence and uniqueness of a Root barrier if [X] = oo. We prepare the proof of Theorem 3 with 
a lemma. 

Lemma 10. Let assumption 1 hold and that 3N 6 N s.t. v\{ m -N,m+N\ = Y^i=iPi°~xi for some n G N and 

(32) v ((-oo, m - N)) = \i ((-oo, to - N)) and v ((to + N, oo)) = p, ((m + N, oo)) . 
Then \1Z reg (er, fj,, v) \ = 1 and E [[X] TR ] = § K xh> (dx) . 

Proof. Define Tjv := [m — N,m + N] and I N := K\int {In)- Condition (32) implies (by Proposition 2) that (at 
least) 

u u = u M on I N . 

In fact (32) first implies the equality in R\Ijv which is then extended to Iff using the continuity of u^,u v . We 
can then conclude that I N C N^' v . By assumption Y^7=i Pi = ^ ~ v {{~ °°) m ~ N)) — v ((m + iV, oo)). 

Since we are looking for a Root barrier, it is clear that X should only stop at the support points {xi : 1 < i < n} 
of v in 7jv (else the stopped process would put mass at a point that is not in the support of v\i N ). Take N^' v 
and define the two indices sets I — {1, ■ • • , n] \ J and J = {j G {1, • • ■ , n} : Xj G iV M,!y }. For some (unknown) 
b = (bi, . . . , b n ) G [0, oo)™ with bj = for any j G ^7, we are then looking for a Root barrier Rf, of the form 

(33) Rb=(^J[bi,+oo]x{xi}J\Jtfi'' v . 

It is clear that if I = (implying that 6 = 0) then the embedding is trivial since u u = iff v = /i. 
Every such barrier Rb gives rise to a first hitting time of the time-space process (t,X t ) t>o 

T b :=mf{t>0:(t,X t )eR b }. 

By assumption [X] 00 = 00 which implies r& < 00 a.s. since T b is dominated by the first exit time of X of the 
interval [xi,x„] (which if finite a.s. since [X] = 00 implies lim sup^^ X t = 00 and liminft^ooXt = —00 see 
for example [29, Proposition V.1.8 ]). 

The subset of barriers r„: Consider c G [0, 00)™ with a corresponding Root barrier R c resp. hitting time r c . 
Define 

(34) r„ = {c G [0, oo)" : cj = Vj G J and P [X Tc = x t ] < p t Vi G 1} 

(for c G r„, X rH embeds "less mass" at points Xi, i G I, than required by i/). Note that if c G T v then one 
has [X Tc = Xj] > ^2j e jPj since one never embeds more mass (under ^) than required at any point in 

the interior (N^' v ) c (i.e. when u M ^ w y ), consequently the excess mass must be embedded at the boundary and 
singular points of N^' v , viz. the points Xj with j G J . 

Y v has a minimal element 7 (in Y v ): We argue as in [21]. We claim that T u with the binary relation "<" 
(component wise) is a lower semi- lattice, i.e. (any finite subset of Y v has a minimal element in Y v ) 

(35) if b, c G r„ then 7 = 6 A c :=(&,• A Cf) f= i,... ,„ G T y . 

To see this, fix i G {1, • • • , n} (for i G ^7 the argument is trivial as 6; = Cj) and assume wlog that bi < Ci. Then 
t 7 = T;, on X T ^ = and t 7 < t& otherwise. Since the barrier points 7 are the smallest of the corresponding 
points between b and c, the hitting time of 7 must be smaller or equal than the hitting times of b and c. Thus, 
{X Ty = Xi} C {X Tb = Xi] and IP [X Ty = a;,-] < p%. The statement follows for general i G {1, • ■ ■ , n}. 

It is clear that the space of barriers of the form of (33) has a minimal element (6 = 0) and being T„ a lower 
semi-lattice it is also clear that an infimum, say d, to Y v exists (wrt to the partial order relation <): it may 
happen though that d £ Y v . Let Td be the first hitting time of Rd, we now prove that 6 I„: if d is the 
infimum to Yy then there exists a monotonic decreasing sequence (dk) ken of elements in Y v converging to d. 
The stopping times Td k induced by the elements dk are uniformly (in k) dominated by the almost surely finite 
first exit time of the barrier A/" M,I/ . Using now Lemma 9 we can conclude that Td k — > Td in probability (and 
hence has an a.s. converging subsequence). By the definition of Y u (see (34)) it follows that 



Vi G X lim P 

k— >oo 



< Pl ^v [x Td = x l ]< Pl d g r„ 



where the implication follows by the convergence of Td k — > Td (a.s. and in Probability) and the continuity of the 
paths of X. 

The minimum ofY v embeds v: Set 7 as the minimal element 7 = (71, . . . ,7„) G Y v (wrt to the partial order 
<) of Y v and associate the induced barrier i? 7 and its hitting time r 7 . We now show that r 7 embeds v via X. 
To see this assume X T v. In this case 

3iei s.t. P [x^ = Xl ] < Pl 

but this contradicts the minimality of 7 since we can find 7 G [0, oo) n by setting 7 fc = 7^ for k ^ i and j i =7^ — 6 
where e > is chosen small enough s.t. P [X T - = Xi\ < pi. Then 7 G T 7 since i? 7 D i? 7 with the only difference 
between the two barriers being the point 7^ — e instead of 7.;. That is, we have increased the probability paths of 
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(t, X t ) hitting the line [7^, +00] x {xi} and therefore decreased the probability that the other lines [7^, 00] x {xj}, 
j 7^ i, are hit, that is P [X T - = Xj] < P [X Ty = xA < pj for j G X \ {i}. Hence, 7 G T v which contradicts the 
minimality of 7. 

The equality E [[X[ T 1 = Jj^x 2 !/ (da;). Write t 7 for the hitting time of i? 7 , 7 being the minimal element of T v . 
Then r 7 is bounded by the first exit time of X from [x\, x n ] which is a.s. finite (see the arguments above before 
we introduced T^). This also implies that t 7 is uniformly integrable. Since r 7 embeds v, it holds E [-X7V.J = 
/ xv (dx) = m. Using that r 7 is uniformly integrable along with the measure v having finite second moment 

and X being a (continuous square-integrable) local martingale one sees E [X] T = E X 2 ^ = J x 2 v (dx) < 00. 

Uniqueness of the barrier (applying Theorem 2). Assume R a , Rb G 1Z (a, fi, v) both of the type (33). Assume 
wlog cti < bi for some i. This implies P [X Ta = Xi] > P [Xt,, = x%] which contradicts that both embed v, 
i.e. v ({an}) = P [X Ta = Xi ] = P [X Tb = Xi \. □ 

The main theorem is proved via a passage to the limit argument (using Lemma 9) where one approximates 
the general distribution v by a sequence of probability measures (fjv) n with finite support in a certain region 
and that converges to v in weakly. We first state a lemma on the existence of the approximating sequence where 
each measure satisfies conditions on the line of that in Lemma 10. 

Lemma 11. Let m € R, /U, v G M 2 n verifying fj, < cx v. Then there exists a sequence of measures (fjv)jveN C 
M 2 m s.t. 

(1) vn converges weakly to v as N — > 00, 

(2) VN\[m-N,m+N]= Z)i=i Pf S x? f or some ( x i ) i=1 C [m - A, m + N], 

(3) l / Ar|]R\[m-JV,m+iV] = H\m.\[m-N,m+N] an d f- <ci V N <cx V . 

Proof. The existence of such a sequence is easy as it follows a construction of the potential functions u UN similar 
to Dubin's construction (see 8 p332-333 in [27]) then via point (5) of Proposition 2 it follows that i>n => v. Fix 
N and have in mind the potential picture of /z, v. Now apply a Chacon- Wash like potential construction but 
starting from instead of us using, A N lines that includes forcefully the points (m — N,u^ (m — N)) and 
(m + A, (m + A)). The obtained function is u VN . 

By construction u VN is a true potential function (simply follow the arguments in [9]), it coincides with in 
1 C N and u„ < u VN < in In . Moreover, the measure i>n induced by u„ N (see Proposition 1) is of finite support 
in Jjv hence i/jy € M-m as well as [i < cx i/jv < C ;c ^- It is also clear that as u VN — >n u v and hence v?q => v (see 
Proposition 2). □ 

We can now prove Theorem 3. In the case 1Z (1, <5o, v) (see [30]) one simply approximates 11 by a sequence of 
finitely supported measures (p- n ) n but for our more general case this is not possible (unless one approximates 
as well the initial measure which makes the proofs much more involved). 

Proof of Theorem 3. For A £ N define ijv := [m — A, m + A] and I N := R \ (In)°- Take a sequence (fjv)jveN 
of probability measures such that i/jv v and where the measures (^jv)jv are as m Lemma 11. 

It is clear that each i>n fulfills the assumptions of Lemma 10 and hence lZ reg (a, fi, vn) = {Rn} and we denote 
with tjv the corresponding hitting time (minimal and uniformly integrable s.th. X TN ~ i/jv and E[[X] TJV ] = 
E [X 2 N ] = J x 2 vn (dx) < J x 2 v (dx) < 00). The compactness of (1Z, p) guarantees the convergence of (Rn) n C 
1Z along some subsequence to some R G 7Z. The problem setup and the embedding properties yield 

sup^^E [[^] TJV ] < K < 00 and [X] «, = 00. 

The continuity of the quadratic variation process implies supjY gN E [rjv] < A' < 00 for some A > (otherwise 
there would exist a subsequence (tn) n with E [rjy] — > 00 and therefore E [[X] Tff ] = J x 2 vn (dx) > K). By the 

second part of Lemma 9 we conclude that E [tr] < K < 00 (as well as E [[X] TR ] < K < 00) and that tn tr 
in probability. 

It remains to show that X TR ~ ^. Now VA G N and V/ G C\ 

f{x) v N (dx)=E[f (X TN )} 

and by the above there exists a subsequence (tn) n s.t. tn Th a-s.. Dominated convergence together 

with continuity of the sample paths of X (and /) along with v show that in the limit V/ G C& 

J R / (x) v (dx) = E [/ (X TR )\ or equivalently that X TR ~ 1^ and hence i? G 7?. re£ , (a, /i, v). The uniqueness of i? 
follows from Theorem 2. □ 



For a more graphical construction with detailed explanation, of which Dubin's solution is a particular case, we point the reader 
to the Chacon- Walsh embedding [9]j see alternatively p337-339 in [27]. 
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7.2. Existence and uniqueness of a Root barrier if X is a GBM. If a(t,x) = x, the process X is a 
GBM which implies [X] < oo and [In X] = oo (this follows for example by Dambis/Dubins-Schwarz), hence 
the results of the previous section do not apply immediately to this important case. Nevertheless, some simple 
modifications allow for a similar argumentation. We prepare the proof by noting that Lemma 9 also applies to 
Geometric Brownian motion. 

Lemma 12. Let a (x) = x, then the conclusions of Lemma 9 remain valid. 

Proof. The proof is a direct reformulation of the proofs of Lemma 2.4 of [30] and Lemma 1 of [2(i] which were 
stated with the Brownian motion being the underlying time-space process. □ 

We can now prove Theorem 4. 

Proof of Theorem 4- The proof follows the same method as the proof of Theorem 3. Assume m = 1. It is clear 
that Lemma 10 holds for this setting (with the proper adaptation of the proof). We point from Assumption 2 
that the support of fj. is bounded from below by a number e. This is crucial as it implies that the hitting time 
to any Root barrier (in this setting) will be dominated by the first exit time of the GBM to the interval [e, oo). 

Define a sequence of approximating probability measures (vn) jvgn such that i/jy => v and define the sequence 
in the same fashion as in the proof of Theorem 3, s.th. Vn is finitely supported in [e, iV] and outside it it 
coincides with fj,. Since each i/jv is finitely supported in [e, N] (outside it we have u M = u VN ) Lemma 10 yields 
the existence of a barrier Rn and its (uniformly integrable) corresponding first hitting time rjv s.th. X TN ~ v n 
and E [[X] TJV ] = J x 2 v N (dx) < K < oo with K > 0. 

From compactness of the barrier space (1Z, r), we can extract a converging subsequence of (Rn) nen with 
limit, say R, and let Tr be its corresponding first hitting time (we abuse notation and still denote by (Rn) n 
the converging subsequence). Now remark that all the stopping times Tjv are dominated by the stopping time 
T e the first time of (X t ) t>o to the interval [e, oo) (see the observations above). Since r e is a.s. finite then so 
is the sequence T/v uniformly dominated and hence Lemma 12 can be applied. The rest of the passage to the 
limit argument follows like the proof of Theorem 3. □ 

Remark 17. If supp (y) C [0, oo) (with v ({0}) = 0) then above argument does not ensure the uniform bound- 
edness of (Tjy) via the uniform boundedness of the L 2 -norms of [X] Tn and this is the reason for Assumption 2 
(i.e. supp (v) C [e, oo) with e > 0). 

7.2.1. Optimality of Root's solution for Geometric Brownian motion. For the general case of embedding a 
probability distribution via an ltd process of the form (6) where a satisfies Assumption 1 with the additional 
assumption of time-homogeneity along with Assumption 2 for /i, v the optimality result has already been stated 
and proved in Rost [34] and Cox- Wang [11]. The particular case of GBM follows also from [11, Theorem 

5.3, Remark 5.5]. We recall their result: in its due generality, Root's construction solves the following general 
minimization problem which is equivalent to Rost's [34] "minimal residual expectation" optimality criterion, 
namely: 



Minimize: E [F (r)] 
subject to: X T ~ v 

t a UI stopping time. 

Here we work under the assumptions of 9 Theorem 4 and F is a given convex, increasing function with bounded 
right derivative / and satisfying F (0) = 0. Under the assumptions of Theorem 4 there exists a Root Barrier R 
and its corresponding stopping time tr which perform the embedding. Define the function 

(36) M(t,x):=E^[f(T R )}, (t, x) e [0, +co) x (0, +co) 

as well as 

" y M (0, z) 



(37) 



Z{x) := 2 



-dzdy x e (0, +oo) . 



Theorem 8 (Cox- Wang [11], Optimality of the Root embedding). Let the assumptions of Theorem 4 hold. Let 
R £ lZ(cr, fj,,!;) and F be any given convex, increasing function satisfying F (0) = and with (bounded) right 
derivative f (x) := F + (x). Assume further that (36) is locally bounded on [0,oo) x (0, +oo) and that for any 
T > the function ( 37) satisfies 

rT 



(38) 



E 



\Z' (X s )a{X s )\ 2 ds 



< oo and E [Z (X )] < 



We point out that the minimization problem here is formulated without the minimality condition, see Definition 6. This is not 
a problem, as under Assumption 2 the state space of X will be in [e, oo), hence a > 0. In this case it is possible to show that the 
uniform integrability of the stopping time is equivalent to the minimality (in the sense of Definition 2). See e.g. Section 8 in [27] 
or Section 3.4 in [21]. 
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Then for any stopping time r for which X T ~ v it holds that 

E [F (tb)] < E [F (r)] . 

Proof. See Theorem 5.3 and Remark 5.5 in [11]. For the particular setting of GBM with initial condition 
Xq ~ Si both (36) and (38) are easily satisfied, see [37, Theorem 4.3.4]. □ 

8. Appendix: comparison theorem 

Comparison theorems for obstacle problems can be found in the literature, see [15] or [23, 22]. However, due 
to the unboundedness of the coefficients they do not cover immediately our setup. For the convenience of the 
reader we provide a complete proof by adapting some results found in [23, 22, 13]. Further, it shows the Holder 
regularity in space of viscosity solutions (which is the first step to get convergence rates for numerical schemes). 

Theorem 9. Let h € C (R, K) be of linear growth, i.e. 3c > such that 

\h(x)\ < c(l + \x\) for all x G R 

and a G C([0,T] x R,R) Lipschitz in space, uniformly in time (swp t \o~ (t, -)\ Lip < oo). Define 

a 2 (t,x) 



Fobs (t, x, r, a, p, M) = min yu (t, x) — h (x) , a 

Let u G USC ([0,T] x R, R) be a viscosity subsolution and v G LSC ([0, T] x R, R) a viscosity supersolution of 
the PDF 

F obs (t, x, d t u, Du, D 2 u) < < Fobs (t, x, d t v, Dv, D 2 v) on (0, T) x R 

Further assume that V(i, x) G [0, T) x R, u(t,x) ,—v(t,x) < C (1 + |x|) for some constant C > and that 
u (0, .) and v (0, .) are 5-Hdlder continuous. Then there exists a constant c > s.t. V (t, x, y) G [0, T) x R x R 

u (t, x) — v (t, y) < sup (uq (z) — vq (z)) + c inf j a~ + a \x — y\ 2 \ . 

Direct consequences of this estimate are 

(1) uq < vq implies u < v on [0,T) x R, 

(2) ifu is also a supersolution (viz. u is a viscosity solution) then u is S-Holder continuous in space uniformly 
in time on [0,T) ; i.e. 

sup \u(t,.)\ c s m < oo. 
te[o,T) 

Proof. Wlog we can replace the parabolic part in F with ^ — a 2 Aw — w (by replacing u resp. v with e~ t u 
resp. e^v). Further we can assume that Ve > 0, u is a subsolution of 

(39) Fobs (t, x, d t u, Du, D 2 u) < - ' 

limu (t, x) = —oo uniformly on O 
(by replacing u with u — jr-r)- Define for a > 0,e > 

i>(t,%, y) = u - v - (j)(t,x,y) with (t, x,y) = e a\x-y\ +el\x\ + \y\ J 

and 

m' e = sup tp (0, x, y) and m Qi£ = sup ijj (t,x,y) — m a e . 

Ixl [0,T]xIxl 

The growth assumptions on u and v together with (39) guarantee for every a > 0, e > the existence of a triple 
(i,x,y) G [0, T) x R x R s.t. 

m a ,e +m° a e = ij)(t,x,y) . 
The proof strategy is classic: the above implies that Va > 0, e > and V (t, x, y) 

(40) u (t, x)-v (t, y) < m Q , e + m° >e + e xt a \x - y\ 2 + e (|x| 2 + |y| 2 ) . 

Using the Holder continuity of uq and we immediately get an upper bound for m° e 

m° aj£ < "o (x) - v (x) + \v \ s \x - yf - a \x - y\ 2 

. j_ 

< \Uo - v \ + ca 2 ~" 

and below we use the parabolic theorem of sums to show that 

(41) m a ^ e < Ca~^ + ke + uj a (e) 
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where ui a (.) is a modulus of continuity for every a > 0. Plugging these two estimates into (40) gives 

u(t,x)-v(t,y) < \u - v \ + {c + C)a-^ + e xt a\x - y\ 2 + e (\x\ 2 + \y\ 2 ^j . 
Now letting e — > and subsequently optimizing over a yields the key estimate 

u(t,x) - v (t,y) < \u Q - v \ + M ((c + C) + e xt a\x - y\ 2 ) . 
Applying it with x = y gives point (1) of our statement. Applying it with a viscosity solution u = v gives 

u (t, x)-u it, y) < inf ( (c + C) + e xt a \x - y\ 2 ) 

and the estimate inf Q> o la" 5 ^ 7 + or 2 j = cr s yields the <5-Holder regularity. 

It remains to show (41). Below we assume m Q:C > and derive the upper bound (41) (which then also holds 
if m a ^ < 0). Note that m Q e > implies i > 0. The parabolic Theorem of sums [12, Theorem 8.3] shows 
existence of 

(a, D x ip (t, x, y) , X) € P 2 ' + u (i, x) and (b, D y ip (i, x, y) ,Y) € V 2 '~v (i, x) 

such that 

(42) a - b = j, (t, x, y) and ( * -y) ^ keUa ( -1 I ) + ke ( \ 

Since u is a subsolution resp. v is a supersolution 

-X -u (t,x) ,u (t,x) - h(x) ] < 



a 2 (t,x) 



min U> - ° ^ ^ Y - v (t, x) , v (i, x) - h (x)j > 

and subtracting the second inequality from the first leads to 

min ^a-b- - X + - ^' ^ Y - u (t, x) + v (t, y) , u (t, x) - v (t, y) -h(x) + h (y)^ < 0. 

First assume the second term in the min is less than or equal to 0. This gives 

u (i,x) - v (i,y) <h(x)-h (y) < \h\ s \x - y\ S 

hence we get the estimate 

(43) m Qje < \h\ s \x - yf . 

Now assume the first term in the min is less than or equal to 0. This gives 

i, (t, x,y)+u (i, x) - v (t, y) < ^S^1 X - ^My 
hence from the definition of ip resp. vn a ^ it then follows that 

\e xt a \x - y\ 2 + m a , e < ^ (a 2 (t, x) X - a 2 (t, y) Y) . 

Estimate the rhs by multiplying the matrix inequality (42) from the left respectively right with the vector 
(cr (i, x) , a (i, y)) resp. (a (i, x) , a (i, y)) to get 

Xe xt a \x - y\ 2 + m a ^ < ke xt a \a (i 7 x) - a (i,y)\ 2 + ke (a 2 (i,x) + a 2 (i,y)) 

(44) < fce A *a|cr|j \x - yf + ke (a 2 (i,x) + a 2 (t, y)) . 
By adding (43) and (44) together and choosing A = lc^fc + 1 we finally get 

m a , £ < \h\ 8 \x - yf - e xt a \x - y\ 2 + ke (a 2 (t, x) + a 2 (t, y)) . 

We estimate the sum of the first two terms on the rhs by using that sup r>0 (r s + jr 2 ^j < ca _5 ^7 and the last 
term using linear growth of a 2 to arrive at 

m Qye < Ca~^ + ke (l + \x\ 2 + \y\ 2 

By lemma (13) we can replace e ^|i| 2 + \y\ 2 ^j by a modulus uj a (e) (i.e. for every a > 0, uj a £ C ([0, oo) , [0, oo)), 
w Q (0) = and u> a is non-decreasing), i.e. 



m a ,e < Ca + ke + ui a (e) . 
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Hence we have shown that Vq > 

limsupm QC = CaT 2 ~ s , 

e 

and 

u (t, x)-v (t, y) < m Q , e + + e At | \x - y\ 2 + e (\x\ 2 + |y| 2 ) . 

□ 

Lemma 13. Let f e USC ([0, T] x M x R) and bounded from above. Set 

m := sup f(t,x,y) 

[0J]xlxl 

m e := sup /(t, x, y) - e (|x| 2 + \y\ 2 ) for e > 0. 

[0,T]xRxR v 7 

Denote with (t e ,x e ,y e ) points where the sup is attained. Then 

(1) lim e ^o m e = sup [0 T]xRxR / (t, x, y), 

(2) e (|i e | 2 + |y e | 2 ) _> e o. 

Proof. By definition of a supremum there exists for every i)>0a triple (t v , x v ,yr]) G [0,T] xlxK such that 
/ (tf/, Xfj, y n ) > m — r/. Fix ?/ > and take e' small enough s.t. e' | 2 + ly^ 2 ^ < ?7- Then Ve G [0, e'\ we have 

m>m e >f (t n ,x n ,y v ) - e [\x n \ 2 + \y n \ 2 ^j > f {t n ,x Vl y n ) -r)>m- 2?/. 

Since r\ can be arbitrary small and e i— > m e is non-increasing the first claim follows. From the above estimate 
and the boundedness of / from above also show that 

k e = e(\x t \ 2 + \y t \ 2 ) 

is bounded. Hence there exists a subsequence of (k e ) e>0 which we denote with slight abuse of notation again as 
(fc e ) e>0 which converges to some limit denoted k{> 0). Now 

f(i € ,x e ,y e ) -k e <m~k e 

and from the first part we can send e to along the subsequence and see that m — k < m, hence k = 0. Since 
we have shown that every subsequence (k e ) converges to the second statement follows. □ 
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